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ABSTRACT 

Based on a homogeneous set of X-ray, infrared and ultraviolet observations from 
Chandra, Spitzer, GALEX and 2MASS archives, we study populations of high-mass 
X-ray binaries (HMXBs) in a sample of 29 nearby star-forming galaxies and their 
relation with the star formation rate (SFR). In agreement with previous results, we 
find that HMXBs are a good tracer of the recent star formation activity in the host 
galaxy and their collective luminosity and number scale with the SFR, in particular, 
Lx ~ 2.6- 10"^^ X SFR. However, the scaling relations still bear a rather large dispersion 
of rms ~ 0.4 dex, which we believe is of a physical origin. 

We present the catalog of 1057 X-ray sources detected within the D25 ellipse for 
galaxies of our sample and construct the average X-ray luminosity function (XLF) 
of HMXBs with substantially improved statistical accuracy and better control of sys- 
tematic effects than achieved in previous studies. The XLF follows a power law with 
slope of 1.6 in the log(Lx) ^ 35 — 40 luminosity range with a moderately significant 
evidence for a break or cut-off at Lx ~ lO^'^erg/s. As before, we did not find any 
features at the Eddington limit for a neutron star or a stellar mass black hole. 

We discuss implications of our results for the theory of binary evolution. In par- 
ticular we estimate the fraction of compact objects that once upon their lifetime 
experienced an X-ray active phase powered by accretion from a high mass companion 
and obtain a rather large number, fx ~ 0.2 x (0.1 Myr/rx) (tx is the life time of the 
X-ray active phase). This is 4 orders of magnitude more frequent than in LMXBs. 
We also derive constrains on the mass distribution of the secondary star in HMXBs. 

Key words: X-rays: binaries - X-rays: galaxies - galaxies: starburst - galaxies: spirals 
- galaxies: irregulars - stars: formation. 



1 INTRODUCTION 

The recent star formation rate of a galaxy is primarily mea- 
sured from the light emitted by massive stars: given their 
short lifetimes, their collective luminosity is, to the first ap- 
proximation, proportional to the rate at which they are cur- 
rently being formed. The scale factors relating luminosities 
at different wavelengths and SFR are usually de rived from 
stellar population synthesis models (see review bv lKennicuttI 
[l998. ). Most of the conventional SFR indicators are known 
to be subject to various uncertainties, one of which is caused 
by the interstellar extinction. 

An alternative method of SFR determi- 
nation is based on the X -ray emission of a 



galaxy (eg. 
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iRanalli. Comastri. fc Settill2003l ). Indeed, X-ray output of 
a normal gas-poor galaxy is dominated by the collective 
e mission of its X-r ay binary (XRB) population (see review 
by lFabbiandl200^ ) . As the evolutionary time scale of high- 
mass X-ray binaries does not excee d ~ few x 10^ yrs, their 
X-ray emission is nearly prompt IShtykovskiy fc Gilfanov! 
I2OO7!). That makes them a potentially good tracer of the 
recent star formation activity in the host galaxy. However, 
in distant galaxies, emission from HMXBs can not be 
separated from that of low-mass X-ray binaries (LMXBs) 
having much longer evolutionary time scales and therefore 
unrelated to recent star-formation or from the contribution 
of the supermassive black hole. For this reason the X-ray 
based method is suitable for young normal galaxies with 
sufficiently large value of the specific, per unit stellar 
mass, SFR, the X-ray output of which is dominated by 
HMXBs. Due to large penetrating power of X-rays, the 
X-ray based method is much less affected by absorption 
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than conventional indicators. This makes it an important 
tool for cross-calibration of different SFR proxies and for 
star-formation diagnostics in galaxies, especially in the 
distant ones. 

Although th e theoretical idea was form ulated 
in the 70-ies (jSunvaev. Tinslev. fc Meier! Il978l l and 
the first supporting ev i dence was obtained in 90-ies 
ijPavid. Jones, fc Formaiil 1 19921 ). the full observational 
justification and calibration of the method required out- 
standing angular resolution and could be achieved only 
with advent of Chandra. During the decade following its 
launch several authors have discussed the X-ray emis- 
sion of compact sou rces and its relation to the SFR 
of the host galaxy (|Grimm. Gi fanov. fc Sunvaev 20031 : 



RanaUi. Comastri. fc Settil l2003l : IPersic fc Rephaelil l2007l : 
Lehmer et al.l |2010| . and several others). All these works 
have established a general agreement regarding the corre- 
lation of the X-ray luminosity with SFR, however there 
was some discrepancy between the scale factors derived 
by different authors. In addition, a notable dispersion in 
the relation was observed. One of the limiting factors of 
previous studies was the heterogeneous nature of the X-ray 
and SFR data. With the large numbers of galaxies now 
available in Chandra, Spitzer and GALEX public data 
archives it has become possible to construct large samples 
of galaxies for which homogeneous sets of multi-wavelength 
measurements are available. This motivated us to revisit 
the problem of the Lx — SFR relation in order to obtain its 
more accurate calibration and understand the origin of the 
dispersion seen in previous studies. 

Study of the populations of HMXBs and its depen- 
dence on the metallicity and other properties of the host 
galaxy is of considerable interest on its own. The specific 
luminosity and numbers of sources as well as their lumi- 
nosity distribution can be used for verification and calibra- 
tion of_th£_bnwj;jevoto synthesis mod- 
els, [g rimm. Gilfanov. fc SunvaevI (|2003l ) showed that in a 
wide range of SFR, the luminosity distribution of HMXBs 
in a galaxy can be approximately described by a univer- 
sal luminosity function, which is a power law with slope of 
~ 1.6, whose normalization is proportional to the SFR. They 
searched for features in the averaged XLF of nearby galax- 
ies, corresponding to the Eddington luminosities of neutron 
stars and black holes, but did not find any. Among other 
reasons, of a more physical nature, they also considered the 
possibility of various systematic effects (e.g. distance uncer- 
tainties) which could smear and dilute the expected features. 
A larger sample of galaxies with well controlled systematic 
uncertainties may allow a substantial progress in this direc- 
tion. 

The present paper is the first of a series of publications 
investigating X-ray emission from star-forming galaxies. It 
will focus on the study of the population of HMXBs and its 
relation to the star formation rate of the host galaxy. The 
structure of the paper is as follows. In Section 2 we present 
the selection criteria and the resulting sample of galaxies. 
The X-ray data analysis methods are described In Section 
3. In Section 4 we study the spatial distribution of com- 
pact sources in our sample galaxies in order to minimize 
the contamination by LMXBs, cosmic X-ray background 
(CXB) and active galactic nuclei (AGN). The analysis of 
far-infrared, near-infrared and ultraviolet data is unfold in 



Section 5. In Section 6 we present the star formation rate 
estimators that we use in the present work. In Section 7 
we study luminosity distributions of HMXBs in individual 
galaxies and their average XLF. In Section 8 we investigate 
the behavior of the collective luminosity of HMXBs and its 
dependence on the star-formation rate. In Section 9 and 10 
we discuss astrophysical implications of the obtained results 
and various aspects of their application. In Section 11 we 
summarize our findings. 



2 THE SAMPLE 

2.1 Primary sample of resolved galaxies 

We started to build the sample of star-forming galaxies from 
the entire "Normal Galaxies" section of the Chandra Data 
Archive. It includes over 300 galaxies of all morphological 
types. We based our selection on the following four criteria. 

i) Hubble type: we selected only late-type galaxies, i.e. 
those having nu merical index of stage along the Hubble se- 
quence T > (|de Vaucouleurs, de Vaucouleurs. fc CorwirJ 
ll976i). According to this definition, SO galaxies were ex- 
cluded since the origin of their far-infrared and ultraviolet 
emission is not associated with recent star formation. This 
first cut reduced the sample to 172 spiral and irregular galax- 
ies. 

a) Specific SFR: as star-forming galaxies contain also 
LMXBs, unrelated to the recent star-formation activ- 
ity, the contribution of the latter will contaminate the 
Lx— SFR relation. Although separation of the high- and 
low- mass XRBs in external galaxies is virtually impossi- 
ble, the LMXBs contribution can be controlled in a sta- 
tistical way. To this end we use the fact that their popu- 
lation i s proport i onal to the stellar mass of the host 
galaxy (iGilfanovl 120041'). whereas H MXBs scale with SFR 
(jGrimm. Gilfanov. fc Sunvae 3[2Q0i), therefore the specific 
SFR (SFR/M*) can serve as a measure of their relative con- 
tributions to the X-ray luminosity of the galaxy. We esti- 
mated the stell ar mass of the galaxies as it is described in 
I Gilfanov! |2004l ). using the T^s-band (2.16 jxm) magnitudes 
from The Two Micron All Sky Survey (2MASS). For the pur- 
pose of the sample selection we used IRAS data to estimated 
SFRs. Based on the fiuxes at 60 |i.m and 100 \un taken from 
IRAS catalogues, we first computed the 42.5 — 122.5 [ im flu x 
(FIR) following iHelou. Soifer. fc Rowan- RobinsonI ljl985[ ). 
Adopting a standard correction factor of 1.7 we converted 
the latter into a 8 — 1000 yi m flu x and computed SFR us- 
ing the calibration of iBelll (|2003l ) (eq. 4) . Taking into ac- 
count the shapes of luminosity distributions for high- and 
low-mass X-ray binaries and typical point sources detection 
sensitivities of Chandra data, we estimated that a threshold 
of SFR/M* > 1 X 10"^° yr~^ would be a reasonable com- 
promise, sufficient for our purpose. After this selection the 
number of galaxies reduced to 88. 

in) Exposure time: we used only those observations in 
which the galaxy is observed sufficiently close to the aim- 



point, with sufficiently long exposure-time (to 



15 ks). 



In case multiple observations were available in the Chan- 
dra archive for a given galaxy, we combined them in order 



Table 1. Primary sample: galaxies resolved by Chandra. 
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Radius of the elliptical galactic region that was used for the analysis (see Sect l4.3| l. 
Star formation rate from Spitzer and GALEX data (see Sect.[6|. 
^ Stellar mass from 2MASS data (see Sect.O. 

f Sensitivity limit in decimal logarithmic units, defined as the luminosity level corresponding to a value of the incompleteness function K{L) = 0.6. 

9 Number of XRBs with luminosity above the sensitivity limit, detected within the galactic region where the HMXBs dominate the compact X-ray source population (Sect. |4]l. 
^ Base-10 logarithm of the total X-ray luminosity of galaxies in the 0.5 — 8keV band, estimated by using eq. I I21II . 
* Predicted number CXB sources within the selected region (see Sect.[4]l, with luminosity above the sensitivity limit. 
^ Predicted contribution of CXB sources to the total X-ray luminosity of galaxies. 

^ for The Antennae we have analyzed a region that includes the two galaxies in the merging system (NGC 4038 and NGC 4039). It is centered at aj2ooo = 180.472 deg and 
'5j2000 = —18.881 deg and has axis ratio=1.4 and PA=78 deg. 
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to improve the sensitivity limit, I/um- After applying this 
selection criterion, 36 galaxies are left in our sample. 

iv) Distance: for the galaxies passed through the above 
criteria, a number of redshift independent distances, deter- 
mined with different methods, were c ollecte d from NEeQ. 
Following the results of iJacobv et al] (|l992l ). we chose the 
most accurate among them. For our primary sample we 
chose a distance threshold of 40 Mpc, in order to spatially 
resolve the X-ray compact sources with Chandra. This also 
allows us to identify and separate the central AGN from the 
rest of detected point sources. After this last selection, 28 
galaxies are left in our sample. 

v) Sensitivity The galaxy list produced in the result of 
this selection had sensitivities in the range of log(Lx) ~ 
35.7 — 38.6 with the median value of log(Lx) ~ 37.1. As we 
perform an accurate incompleteness correction, no further 
filtering based on the point source detection sensitivity was 
performed. 

For the purpose of binary population studies, this rigor- 
ously defined sample was complemented by two more galax- 
ies. In order to improve statistics of the high luminosity 
end of the XLF we added the Cartwheel galaxy. It has the 
largest SFR value among our resolved galaxies and although 
is more distant (D~ 123 Mpc), the brightest sources, with 
log(I/x) 39, do not suffer from significant sources confu- 
sion in Chandra data. To expand the low lu minosity end we 
used the ASCA data of SMC observations (jYokogawa et al.l 
l2003f ). The SMC data were used only for the XLF construc- 
tion. 

The so defined primary sample is presented in Table [T] 
It includes 29 galaxies, with SFR values ranging from ~ 0.1 
to - 20 Mq yr-\ 



2.2 High-SFR sample 

In order to explore the Lx — SFR relation in the high SFR 
regime, we relaxed the distance constraint and constructed a 
secondary sample from the galaxies that passed the first and 
second criteria, and are characterized by high star-formation 
rates, SFR > 15 Mq /yi. In total, we found 8 galaxies with 
SFR ranging from 15 to 290 Mq yr~^ (all, naturally, located 
outside Z) > 40 Mpc). 

These galaxies were treated as spatially unresolved in 
our analysis (see Sect. 17.11 for details and caveats). Five 
of them are classified as ultra-luminous infrared galaxies 
(ULIRGs) and two as luminous infrared galaxies (LIRGs). 
The remaining object (NGC 4676 or The Mice), is a merging 
system. 

For these galaxies we verified that the AGN emis- 
sion does not compromise the IR-b ase SFR det e rmina - 
tion. To this end, we used results of Farrah et al.l ll2003h ; 
iNardini et all l|2009l '): lYuan. Kewlev. fc Sanders! (|2010l ') who 
used various methods (IR and optical photometry, 5 — 8 ^ 
Spitzer spectroscopy and SDSS data) to constrain the con- 
tribution of the AGN to the bolometric luminosity. Their 
results are summarized in the last two columns of the Ta- 
ble [2] In all but one galaxy the AGN contribution to the 
bolometric luminosity is constrained to < 30%. The remain- 
ing galaxy, IRAS 13362-^483 (NGC5256), is classified as a 



2009). However, the [NV]/[NII] 



Seyfert 2 jNardini et all 

line ratios ( Petric et al.l 20111 ) and absence of the broad 
line region ( Wu et al.l I2OIII ') in this galaxy suggest that it 
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is starburst-dominated. We therefore decided to keep it in 
our sample. 

Finally, we also included 8 normal galaxies from 
the Chandra Deep Fiel d North (HDFN) l|Brandt et all 
2001) an d Lynx field JStern et all I2OO2I ) compiled by 
Grimm, Gilfanov. fc Sunvae vl l|2003l ). Their selection is 
described in the above mentioned paper; their K- 
corrected X-ray and radio l u minos ities were taken from 
Gilf anov. Grimm, fc SunvaevI (|2004l ). Thus, our high- SFR 
sample included 16 galaxies, listed in Table [5] and O 



2.3 Observer bias 

Obviously, the selection of galaxies based on the content of 
the Chandra archive is subject to the observer bias - we 
often tend to choose galaxies for our proposals which are 
known to be bright in X-rays. This can shift the scale factor 
in the obtained Lx — SFR relation upwards. It may also 
affect the shape of the luminosity function and that of the 
Lx — SFR relation. As it is further discussed in the Section 
18.31 the behavior of the latter suggests that the observer 
bias does have some effect on our sample. A straightforward 
method to assess its amplitude is to compare the Lx — SFR 
relation obtained here with that based on galaxy samples 
selected from multi- wavelength surveys. This work, although 
important, is beyond the scope of this paper and will be 
addressed in a follow-up publication. 



2.4 A remark regarding the roles of primary and 
secondary samples 

It is important to realize the different roles and the differ- 
ence in the way the two samples are treated. The main goal 
of this paper is to study population of high-mass X-ray bi- 
naries, their luminosity distribution and dependence of their 
number and collective luminosity on the star-formation rate 
and other properties of the host galaxy. This is achieved 
with the help of resolved galaxies of our primary sample. 
Their analysis is presently limited to the emission from com- 
pact sources, the contribution of various unresolved emission 
components is not considered. To properly resolve the pop- 
ulations of compact sources, we considered only galaxies lo- 
cated within the distance of L) < 40 Mpc. This limited the 
volume and, correspondingly the range of SFRs probed by 
the primary sample. Indeed, the highest SFR in our resolved 
galaxies does not exceed ~ 15 — 20 Mq/jt. 

In order to probe the high-SFR regime, we comple- 
mented our primary sample with more distant galaxies. This 
allowed us to extend the dynamical range of SFRs by ~ 
two orders of magnitude. On the other hand, due to larger 
distances to these galaxies, their compact source popula- 
tions can not be studied at the same level of detail, because 
of insufficient point source detection sensitivity and source 
confusion, even in the most nearby of them. Therefore we 
treated them as unresolved (even though the angular ex- 
tent of some of them was as large as ~ 25 — 50") and only 
measured their total X-ray emission. Obviously, X-ray lu- 
minosities of these galaxies can not be directly compared 



HMXBs in star-forming galaxies 5 



Table 2. High-SFR sample: LIRGs and ULIRGs. 
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distance from lSanders et al ] l|2003l l. 
^ AGN contribution to the bolometric luminosity (in per cent). 

Ref eren ces for AGN contribution to th e bolome tric luminosity: ( l lNardini et al.l l l2009tl , (2 lFarrah et al.l 
l|2003l ). C3 lYuan. Kewlev. fc Sanderj l|2010 '). l'4'Bcr nard-SaIas et alj (2009l 'l 



Table 3. High-SFR sample: galaxies from the Hubble Deep Field North and Lynx Field. 



Galaxy 


Redshift 


-fb.5-8kcV 


io.S-SkoV 


i^l.4GHz 


il.4GHz 


SFR 






(10-1^ erg cm^^) 


(10*0 gj.g/g) 


( ^tJy) 


(1030 erg/s) 


(Mq yr-i) 


123634.5-1-621213 


0.456 


0.43 


31.43 


233.00 


1.63 


90.56 


123634.5-1-621241 


1.219 


0.31 


163.86 


230.00 


15.94 


885.56 


123649.7-1-621313 


0.475 


0.15 


10.09 


49.20 


0.37 


20.56 


123651.14-621030 


0.410 


0.30 


14.61 


95.00 


0.52 


28.89 


123653.44-621139 


1.275 


0.22 


128.03 


65.70 


5.18 


287.78 


123708.34-621055 


0.423 


0.18 


9.39 


45.10 


0.23 


12.78 


123716.34-621512 


0.232 


0.18 


2.55 


187.00 


0.27 


15.0 


0084857.7+445608 


0.622 


1.46 


177.76 


320.00 


4.50 


250.0 



Note: See iGrimm. Gilfanov. fc Suny acv ( 2003 ) for a complete d e script ion of the sample; the X-ray and 
radio luminosities are adopted from lGilfanov. Grimm, &: SunvaevI l l2004h . 



with our primary sample, as they include the contribution 
of diffuse emission, which may be significant in star-burst 
galaxies. However these galaxies help to present the broader 
picture and to constrain the properties of bright HMXBs and 
ULXs in the high SFR regime and at intermediate redshifts, 
based on the beh avior of their collective emission, usin g the 
ideas described in lGilfanov. Grimm, fc SunvaevI (|2004 ). The 
quantitative analysis of the full Lx — SFR relation for the 
two samples, including diffuse X-ray component, will be pre- 
sented in the forthcoming publication. 

3 X-RAY ANALYSIS 
3.1 Data preparation 

We analyzed in total 64 Chandra observations, 63 ACIS-S 
and one ACIS-I, listed in TableU) The data preparation was 
done following the standard CIAC0 threads (CIAO version 
3.4; CALDB version 3.4.1) and limiting the energy range 
to 0.5 — 8.0 keV. We performed the point source detection 
in each observation using CIAO wavdetect. The scale pa- 
rameter was changed from the default value. We used the 
-\/2-series from 1.0 to 8.0 to have a wide enough range of 

^ http:/ /cxc. harvard. edu/ciao3.4/index. html 



source sizes to account for the variation in PSF from the 
inner to the outer parts of an observation. In order to avoid 
false detections, the value of the parameter sighthresh was 
kept rather low. We set it as the inverse of the total number 
of pixel in the image. The typical values being in the range 
~ lO"'^ - 10-**. We used maxiter = 10, iterstop = 0.00001 
and bkgsigthresh = 0.0001. We set the parameter eenergy 
= 0.8 (the encircled fraction of source energy used for source 
parameter estimation). 

Readout streaks caused by a bright point source were 
seen in two cases: observation 378 (M82) and 315 (Anten- 
nae). To avoid that such a problem could affect source de- 
tection, we corrected the events using acisreadcorr. 

Multiple observations were present for ~ 30% of the 
sample galaxies; they were combined in order to improve the 
sensitivity and, in case of large objects, to cover the galactic 
area. We first used CIAO reproject_events to reproject 
all the observations into the sky coordinates of the longest 
observation of a given galaxy. The images were then merged 
using merge_all script and the wavdetect task was applied 
again to the combined image. An exposure map weighted 
according to a specific model for the incident spectrum was 
computed for each observation. 

In computing the exposure maps and counts-to-ergs 
conversion factor we assume an absorbed power law spec- 
trum with the photon index F = 2.0. This choice was moti- 
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vated by the value of the centroid of the power law photon 
index distribution of the luminous co mpact X-ray sources 
in star- forming galaxies measured bv ISwartz et al.l l|2004h . 
r = 1.97 ± 0.11. In order to chose the absorption, we pro- 
duced the combined spectrum of all compact sources in 8 
representative galaxies from our sample and determined the 
best fit value of uh ~ 3 x 10^^ cm~^. This value used for all 
galaxies. We note that moderate variations of these parame- 
ters do not result in significant changes of the counts-to-ergs 
conversion coefficient - change of the slope by AF — 0.2 or 
of the absorption by a factor of 1.5 results in less than 10% 
change of the conversion factor. This is insignificant for the 
purpose of our study. 



3.2 Source counts and luminosities 

For X-ray photometry of compact sources we u sed the ap- 
proach and scripts from lVoss fc Gilfanov! (|2007l ). The count 
rate for each detected point source was calculated inside a 
circular region centered on the source central coordinates 
given by wavdetect output. In order to determine the ra- 
dius of the circle, for each observation we extracted the point 
spread function (PSF) using CIAO mkpsf task. In case of 
multiple observations the PSFs were combined using the val- 
ues of the exposure maps as weights. We mapped the PSFs 
into the World Coordinate System (WCS) reference frame 
of the relative point source image using reproject_image 
task. The radius of the circle was determined individually 
for each source so that the encircled PSF energy was 85%. 
For the background region we used an annulus which in- 
ner radius was equal to the radius of the source region and 
the outer radius was 3 times bigger. The corrected source 
counts and errors were then found by the following equa- 
tions (|Voss fc Gilfanov! [2OO7I ): 

C{b-d)d-^ -Q 



S ■ 



2 



abd-^ -13 
al(b~dfd-^ + ol 



(1) 



(2) 



(aM-i - 

where S is the number of net counts from the source, C is 
the number of counts inside the source region and Q is the 
number of counts in the background region, a is the integral 
of the PSF over the source region, /3 is the integral of the 
PSF over the source and background regions, b is the area 
of the source and background regions and d is the area of 
the source region. The latter two quantities were derived by 
integrating the exposure maps rather than computing the 
geometric areas. 

Strictly speaking, the PSF weighting should be based 
on source counts rather than exposure map values. Indeed, 
the compact sources under study may show significant vari- 
ability between observations which may result in the inac- 
curate determination of the shape of the combined PSF, of 
the radius of the 85% region and, consequently, inaccurate 
flux calculation. However, several factors limit the practical 
usefulness of counts-based PSF weighting. Firstly, co-adding 
observations leads to the increase of the sensitivity, there- 
fore a fraction of sources detected in the summed image are 
not detected in individual observations. The fraction of such 
sources is large for the shape of the HMXB luminosity func- 
tion. Secondly, as the source counts are subject to Poissonian 



uncertainty, using them as weights in determining the effec- 
tive PSF will introduce additional statistical uncertainty in 
the source intensity calculation. For faint sources, this ad- 
ditional uncertainty will exceed the inaccuracy of the expo- 
sure map-based procedure. Given the pattern of the Chan- 
dra pointings for the galaxies in our sample, the inaccuracy 
of the exposure-map-based weighting method is relatively 
small. Indeed, the galaxies in our sample were the main 
target of Chandra observations and therefore were similarly 
located in the fov of individual observations. With a few ex- 
ceptions, the analysis regions were located entirely inside a 
single ccd chip. For these reasons, for a given source, there is 
not much variation in the size of the 85% circle between indi- 
vidual observations. To verify this, we selected several bright 
sources located in galaxies with multiple pointings and com- 
pared their intensities obtained using the two psf weighting 
techniques. We found that the difference between the two 
methods did not exceed ~ 5 — 10% (a part of which is due 
to the additional statistical uncertainty introduced by the 
counts-based weighting method). For these reasons, for the 
purpose and parameters of the present work the exposure 
map -based PSF weighting technique is more appropriate. 

The problem of point source crowding may compro- 
mise the total X-ray luminosity estimation of galaxies. We 
identified in each galaxy all the compact sources having 
background regions overlapping neighboring sources. For the 
galaxies affected by crowding problem, we performed a sec- 
ond iteration of the circular photometry procedure defining 
the radius of each source so that it included 90% of encir- 
cled PSF. These regions were then excluded from both the 
individual image and exposure map of one of the overlap- 
ping point sources in order to subtract the source counts and 
correct the source area respectively. We finally used the cor- 
rected image and exposure map to perform the photometry 
procedure defining the radius of the circle for each source 
so that the encircled PSF energy was 85% and obtained the 
corrected C, Q, b and d for applying eqs. [1] and |2l 

The X-ray luminosities of X-ray binaries in galaxies 
from our primary sample were computed as a sum of lu- 
minosities of compact sources, corrected for the incomplete- 
ness and CXB contribution and transformed to a common 
luminosity limit, as described in Section IS] 

The galaxies from the secondary high-SFR sample were 
treated as unresolved objects and their luminosity was mea- 
sured as the background-subtracted luminosity inside the 
D25 ellipse. The counts-to-ergs conversion in this case was 
determined based on the modeling of their X-ray spec- 
tra. For most of the spectra the best-fit was obtained 
with a two component model: a thermal plasma plus a 
power-law corrected for the Galactic absorption. In two 
cases (IRAS 20551-4250 and IRAS 23128-5919) a good 
fit was obtain ed with a thre e-component model, according 
to [Pranccschi ni et al.! (jiooi): a thermal component and a 
"leaky -absorber" continuum, including an absorbed plus 
a non-absorbed power law spectrum of the same photon in- 
dex. To measure the X-ray luminosity of unresolved galaxies 
and the contribution of diffuse gas in resolved galaxies, we 
created source and background spectra and the associated 
ARF and RMF files using specextract script, which is suit- 
able for extended sources. Spectra were then modeled using 
XSPEC V. 12.3.1X. 

The luminosities of HDFN galaxy sample were mea- 
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Table 4. Chandra observations analyzed. 



Galaxy 


Obs-ID 


Galaxy 


Obs-ID 


NGC 0278 


2055 


NGC 


4631 


797 




2056 


NGC 


4676 


2043 


NGC 0520 


2924 


NGC 


5194 (M51A) 


3932 


NGC 1313 


2950 






1622 


NGC 1427A 


3949 


NGC 


5253 


7153 


NGC 1569 


782 






7154 


NGC 1672 


5932 






2032 


NGC 2139 


8196 


NGC 


5457 (MlOl) 


934 


NGC 2403 


4630 






4732 




4628 






5309 




4629 






5322 




2014 






6114 


NGC 3034 {M82) 


378 






4736 




2933 






4731 


NGC 3079 


2038 






5340 


NGC 3310 


2939 






5300 


NGC 3556 


2025 


NGC 


5474 


9546 


NGC 3631 


3951 


NGC 


5775 


2940 


NGC 4038/39 


315 


NGC 


7090 


7252 




3040 






7060 




3041 


NGC 


7541 


7070 




3042 


NGC 


7793 


3954 




3043 


IRAS 17208-0014 


2035 




3044 


IRAS 20551-4250 


2036 




3718 


IRAS 23128-5919 


2037 


NGC 4194 


7071 


IRAS 10565-f2448 


3952 


NGC 4214 


5197 


IRAS 13362-f4831 


2044 




4743 


IRAS 09320-1-6134 


2033 




2030 


IRAS 00344-3349 


8175 


NGC 4490 


4726 


CARTWHEEL 


2019 




4725 


UGC 05720 


9519 




1579 








NGC 4625 


9549 









sured in 0.5 — 8keV band, u si riR th eir count rates, best 
fit slopes from [ Brandt et al.l l|200ll ) and redshifts from 
ICohen et al.l \2Qm . 



4.1 Contribution of tlie central AGN 

In tlie resolved galaxies [D < 40 Mpc), the AGN contri- 
bution can be easily discriminated. Three galaxies from the 
primary sample are known to have a low luminosity AGN: 
NGC 1672, NGC 3079 and M51A ( Vcron-Cctty fc VeronI 
120061 ). Whereas the central sources of M51A, NGC 3079 
were detected and excluded from the further analysis, no 
central X-ray sourc e was found in NGC 1672. However, 
IJenkins et ahl l|201ll ) found in this galaxy an evidence for 
a low-luminosity nucleus, with the 2—10 keV luminosity of 
~4xl0^® erg s"^ Its non-detection in our analysis is likely 
a result of different wavdetect settings, in particular of a 
higher detection threshold used (in fact, visual inspection of 
the image did reveal a hump on top of the diffuse emission in 
the center of the galaxy). In addition to these three galaxies 
we also decided to exclude the compact X-ray source in the 
center of NG C 4039, although its nature may be somewhat 
controversial. iBrandl et al.l l|2009l ) ruled out the presence of 
an AGN in this galaxy, based on the lack o f [A^'eV] emission 
and th e low [OIV]/[NeII] ratio. Moreover. iDesroches fc Hoi 
l|2009D suggested that the X-ray spectrum is indicative of an 
X-ray binary in the very high state, rather than of an AGN. 
However, they pointed out that location of this source in 
the nuclear star cluster makes it a likely LMXB candidate, 
rather than an HMXB. 

For galaxies of the high-SFR sample, the AGN con- 
tribution in the X-ray band was constrained directly from 
the Chandra images. Galaxies from the "unresolved" sam- 
ple of Table [2] are in fact partly resolved by Chandra, the 
smallest one having angular size of D25 ~ 24 arcsec. This 
is enough to verify, that none of them has a dominating 
central point source. The upper limit of 30% on the AGN 
contribution to the IR luminosity (see Sect. [2} ensures that 
the IR based SFR determination is not compromised by the 
AGN. The HDFN and Lynx field galaxies are spiral and ir- 
re gular galaxies, showing no sign s of AG N activity, selected 
bv lCrimm. Gilfanov. fc Sunv acv' (2003') based on the iden- 
tifications from [Richards et al., (^199^). 



4 SPATIAL ANALYSIS 

Whereas irregular galaxies host young stars, the stellar pop- 
ulation of spirals is spread in a broad range of ages. The 
old populations of bulges are totally dominated by LMXBs, 
similar to the ones of elliptical galaxies. An active nucleus 
can be also present. The disk instead hosts X-ray binaries of 
different ages. As galaxies are mostly transparent to X-rays 
above 2 keV, background X-ray sources are also detected in 
each galactic area. Although reliable separation of different 
classes of sources in external galaxies based on their X-ray 
properties is virtually impossible, contributions of LMXBs 
and background AGN can be controlled and subtracted in 
the statistical manner. In order to identify regions of galax- 
ies where the HMXBs dominate the compact X-ray source 
population, we analyzed the spatial distribution of the de- 
tected compact sources in resolved galaxies of the primary 
sample. An example of regions selected based on this method 
is shown in Fig. (2] 



4.2 Contribution of LMXBs 

Although only star-forming galaxies with SFR/Af* ^ 1 x 
10"^" yr~^ have been selected, the contribution of LMXBs 
to the total X-ray luminosity can be further minimized by 
excluding bulge sources, which populations are dominated 
by LMXBs. This procedure was applied to spiral galaxies 
having an inclination angle i < 65 deg where the bulge region 
was identified by inspecting the optical and near-infrared 
images and taking into account the information available 
in literature. The following spiral galaxies were affected: 
NGC 0278 (iibuigo = 0.1 arcmin), NGC 3631 (i^buigc = 0.3 
arcmin), M51 A (i?buigc = 0.4 arcmin), MlOl (i?buigc = 1-5 
arcmin) and NGC 4490 (iibuigc = 0.3 arcmin). For MlOl 
the bulge region was chosen based on results from I Gilfanov! 
(2004). The X-ray point sources detected in such regions 
were then excluded from the analysis. We excluded the same 
regions from the analysis of infrared and ultraviolet data 
(Sect. [5]). The fraction of galaxy removed by this procedure 
ranges between 10% and 30%. 
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4.3 Contribution of CXB sources 

In tlie central regions of galaxies the surface density of 
XRBs is sufficiently high that the contribution of back- 
ground AGNs can be neglected. In the outer parts the sur- 
face density of compact sources becomes comparable to the 
average density of CXB sources and the latter may become 
and important source of contamination. In order to minimize 
this contamination we used the following procedure. 

For each resolved galaxy we defined an elliptical region 
according to the isophotal diameter (D25), isophotal diam- 
eter ratio {R25) and the position an gle of the major axis 
(PA\ taken from the RC3 catalog (|de Vaucouleurs et al.l 
[l993)- We computed the incompleteness function K{L) of 
the part of the Chandra image within this re gion, following 
the method and code from .Voss fc Gilfanovl l|2006l ). We de- 
fined the completeness luminosity Lcomp as the luminosity at 
which K{L) = 0.8 (i.e no more than 20% of point sources are 
missing). For point sources detected with luminosity above 
the completeness limit, we constructed cumulative and dif- 
ferential radial profiles (see Fig. [T] for an example) account- 
ing for the elliptical shape of the region. The latter were 
compared with the CXB level predicted above the same lu- 
minosity threshold . We used the full band log A'' — log S of 
iGeorgakakis et al] l|2008l) and converted it to the 0.5 — 8 keV 
band used for the detection. We checked its consistency with 
the observed density of point sources using differential pro- 
files (left panel in Fig[T]). The radius i?x of the (elliptical) 
galactic region that will be used for the analysis was de- 
fined based on the cumulative profile as the radius within 
which the contribution of CXB sources equals 30% of the to- 
tal number of detected compact sources. The so determined 
radii i?x approximately follow the relation: 

iix(arcmin) = iif ~ 30%^ = 1.2 x (3) 

\ A'tot / \ 2 J 

In four galaxies (NGC 1427A, NGC 4194, NGC 4625 
and UGC 05720) we detected a rather small number (4 — 6) 
of point sources that did not allow us to perform spatial 
analysis. For these galaxies we analyzed the D25 ellipse re- 
gions, as their CXB contamination was found to be in the 
5% - 16% range. 



5 MULTIWAVELENGTH ANALYSIS 

For SFR estimations we used data from Spitzer and GALEX 
archives, the stellar mass was estimated from 2MASS im- 
ages. For the galaxies from the resolved sample we computed 
SFR and mass inside the analysis regions defined in the pre- 
vious section. For the galaxies of high-SFR sample we did 
the measurements in the D25 regions. 

5.1 Far-infrared 

We analyzed MIPS 70 [im and 24 \im Large Field images. In 
particular, we used the "post Basic Calibrated Data" prod- 
ucts provided by the Spitzer Space Telescope Data archival. 
These are images calibrated in MJy/sr, suitable for pho- 
tometric measurements. We did not analyze MIPS 160 \im 

^ http:/ /irsa. ipac.caltech.edu/applications/Spitzer/Spitzer/ 
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Figure 1. An example of the spatial analysis performed in re- 
solved galaxies. Left: surface density of the detected point sources 
(soli d histogram) plotted toge ther with the predicted CXB level 
from lCeoreakakis et al.l (120081) (dotted-dashed line). Right: cumu- 
lative radial profile of detected point sources (solid line) and the 
CXB growth curve (dotted-dashed line). The vertical dashed line 
indicates the radius where the cumulative number of CXB sources 
equals 1/3 of the total number of detected compact sources. It is 
the radius i?x of the elliptical region we choose to analyze. 




Figure 2. Optical (HST) image of MlOl. The crosses indicate the 
positions of detected X-ray point sources. The sources marked in 
green are used for the further analysis. The part of the galaxy 
inside the red circle (X-ray sources marked by red crosses) were 
excluded because their p opulations were demonstrated to be dom- 
inated by LMXBs (see iGilfanml |2004 for details). The region 
outside green circle (X-ray sources marked by cyan crosses) was 
excluded because the contribution of CXB sources becomes large. 
Its size was chosen as described in Sect [4] 



images because for most of the galaxies in our sample only 
"Small Field" images were available. This operating mode 
often covers only a fraction of galaxy. For three galaxies 
(NGC 7541, NGC 2139 and NGC 3631) neither MIPS 70 \im 
nor 24 ^m data were available. For NGC 3034 (M82) Spitzer 
data are affected by saturation effects that limit our ability 
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to extract reliable global f lux densities. In th ese cases we 
used published IRAS data jSanders et al.|[2003l ) to compute 
the IR luminosity. 

After having subtracted the background counts, the to- 
tal net counts were converted from units of MJy/sr into Jy 
using the following conversion factors, which take into ac- 
count the pixel size of MIPS 70 ^.m and 24 \xxa. detectors: 



Croum = 3.76 X 10' 



C24^m = 1.41 X 10' 



(4) 



The monochromatic fluxes (Jy) at 24 \im and 70 |J.m were 
then converted into spectral luminosities (erg/s/Hz). The 
total IR lum inosity (8—1000 \ im) w as estimated using the re- 
lations from lBavouzet et al.l l|2008h . which give Lir in units 
of solar luminosity Lq = 3.839 x 10^^ erg/s.: 



Lir(L0) = 7.90 X (i^L.)?o ; 



Lir(Lq) = 6856 X [vL,) 



24 p.m,rcst 



(5) 
(6) 



We preferentially used equation ( [5ll, as it has a ~ 3 ti mes 
smaller dispersion than eq. ((6)1 (see Bavouzet et al.|[2008l . for 
details) . 



5.2 Near-infrared 

The near-infrared analysis was performed using the Ks 
(2.16 \xm) images and mosaics of 2MASS All-Sky Extended 
Source Catalog] (XSC) and Large Galaxy Atlaffl (LGA). 
We converted the source counts measured on images into 
calibrated magnitudes using the zero point magnitude 
KMAGZP given in the image header: 



mK(mag) = KMAGZP - 2.51og(5') 



(7) 



where S is the integrated, background-subtracted flux in 
data- number units ("DN") measured within the selected 
galactic region. After having obtained magni tudes, we 
estima ted the Tfs-band luminosity and, following iGilfanovl 
l|2004h , the stellar mass. 



5.3 Ultraviolet 

We based the ultraviolet analysis on GALEX NUV and FUV 
background-subtracted intensity maps publicly available in 
the archive of GR4/GR5 Data Releas^fl These are images 
calibrated in units of counts per pixel per second, corrected 
for the relative response, with the background removed. We 
converted the net counts to flux using the conversion factor^ 
between GALEX count rate (ct/s) and flux (erg/cm^ /s/ A): 



Cfuv = 1.40 X 10' 



Cnuv = 2.06 X 10' 



(8) 



The fluxes were then converted into monochromatic lumi- 
nosities (erg/s/A) at 1529 A (FUV) and 2312 A (NUV). For 
three galaxies of our sample (NGC 0278, NGC 4194 and 
NGC 7541) neither NUV nor FUV data were available. Their 
treatment is discussed later. 



* http:/ /irsa. ipac.caItech.edu/applications/2MASS/PubGalPS/ 
^ http: / /irsa. ipac.caltech.edu/applications/2MASS/LGA/ 
^ http: / /galex.stsci.edu/GR4/?page=mastform 
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Figure 3. The SFR — Af* plane. Different types of galaxies are 
plotted with different symbols. The dashed lines correspond to 
constant stellar-mass-to-SFR ratio. 



6 STAR FORMATION RATE 

We used a complex SFR proxy that takes into account 
both the UV light escaping the galaxy and the IR 
emission of the dust heated by young stars LB cU 2003 ; 
Hirashita, Buat. fc Inouell20oj : llglesias-Paramo et al.ll20o3 . 



20061 ). It has the advantage of being free of the model de- 
pendence for the attenuation correction of the UV emission. 
It contains also a correction factor, 77, accounting for the 
fraction of the IR emission due to old stellar population in 
normal star-forming galaxies: 

SFRtot = SFRSuv + (1 - »?)SFRiR (9) 

where SFRg,uv and SFRia are obtained assuming a Salpeter 
IMF from 0.1 to 100 Mq: 

SFRiR(Af0yr-i) = 4.6 x 10-**LiR(erg s"') (10) 



SFR^ 



,(M0yr" 



1.2 X 10" 



^(ergs ^) (11) 



^Nuv.obs is the observed NUV (2312 A) luminosity, i.e. un- 
corrected for dust attenuation, and Lir is the 8— 1000 [ xm lu - 
minosity. Eqs. and (fTTj) bvligTe sias-Par amo et al.l (120061 ) 
are consistent with those of lKennic utt (1998j) within 3% and 
10% respectively. T he second t erm in eq. is very similar 
to the cahbration of lBelll (|2003l ) (their eq. (5)). The quantity 
77 is the cirrus correction, i.e. the fraction of IR luminosity 
due to old stellar population. It depends on the type of the 
galaxy but accurate values for individual galaxies are difli- 
cult to determine , therefore average values are usu ally used. 
We use results of iHirashita. Buat. fc Inoud (|2003h : 



0.4 for normal disk galaxies 
for starbursts 



(12) 



To be consistent with [ Hirashita. Buat."fc^ Inoud l|2003l ), we 



^ http://galexgi.gsfc.nasa.gov/docs/galex/FAQ/counts_background.htnmlsed the atlas of iKinnev et al.l 1 19931 ) to classify the objects 
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Figure 4. Cumulative X-ray luminosity functions of the galaxies 
from primary sample, normalized by their respective SFRs. The 
solid line is the cumulative XLF per unit SFR, given by integra- 
tion of the equation I I18II . 



of our sample as starbursts or normal star-forming galaxies 
and use the appropriate value of rj in computing the SFR. 
Th is de f inition is similar, but not identical to the one used 
by iBelll (|2003l ) (r? ~ 0.32 for galaxies having Lir > IO^Lq 
and ri ~ 0.09 for galaxies having Lir ^ 10^^ Lq). 

For each of the three galaxies with no UV data avail- 
able (see Sect. I5.3|l . we determined the average value of 
the LNUv,obs/LiR ratio for galaxies of similar Hubble type 
and inclination. We obtained values in the range 0.08 — 0.4. 
This ratio was used to estimate I/Nuv.obs from Lir and then 
eq.lO was applied. 

The SFRs for HDFN galaxies were com puted based on 
their 1.4 GHz luminosities and calibration of iBelll l|2003l ): 



SFR(M0 yr" 



5.55 X 10 Li.4GHz(ergs 



(13) 



In Fig. Owe present the two samples of star-forming galaxies 
in the SFR - M, plane. 



7 THE LUMINOSITY FUNCTION OF HMXB 

The Fig. |4] shows cumulative luminosity distributions for 
all galaxies from the primary sample, normalized to their 
respective SFRs. It is apparent that although shapes of the 
distributions are similar, there is a considerable dispersion 
in their normalization. Their amplitude and significance will 
be discussed later in this section. 



7.1 Average XLF of HMXBs 

We construct the average luminosity function combining the 
data for all resolved galaxies. It includes over 700 compact 
sources. An intuitive and straightforward method would be 
to bin the sources into e.g. logarithmically spaced bins and 



normalized the result by the sum of the SFR of all galaxies 
contributing to the given luminosity bins. The latter step 
is required to account for the fact that different galaxies 
have different point source detection sensitivity. However, 
because of the considerable dispersion in the normalization 
(Fig|3)), the so computed luminosity function may have a 
number of artificial steps and features at the luminosities, 
corresponding to the luminosity limits of particular galaxies. 
In order to deal with this problem we used the following 
method. 

Considering that there is much larger dispersion in nor- 
malization than in the shape of the XLF in individual galax- 
ies, we write for the luminosity distribution in the k-th 
galaxy: 



dN 
dL 



= a- SFRfc /(L) 



(14) 



where SFRk is star-formation rate in the k-th galaxy, - 
the XLF normalization and the function /(L) describes the 
XLF shape, assumed to be same in all galaxies, which we 
would like to determine. The number of sources in the j-th 
luminosity bin, AN{Lj), is: 

AN{L,) = J2 AiVfe(L,) = /(L,) AL, ^ SFR^^fe (15) 

where ANk{Lj) is the CXB-subtracted number of compact 
sources in the k-th galaxy in the j-th luminosity bin and 
summation is performed over all galaxies of the sample. The 
/(L) can be determined as: 

1 AiV, 



For a power law luminosity function ,/(L) = L 
are calculated as: 



6 



1-7 Affe(> Lth.fc) -^cxb(> Llh.fc) 



SFRfc 



-L, 



-7+l^ 



(16) 
the 5fc 



(17) 



where Lth.fc is the sensitivity limit for the k-th galaxy, 
Nk{> Lth.fc) is the number of sources detected above this 
sensitivity limit, A'^cxb(> Fth,k) is the predicted number 
of resolved CXB sources above the corresponding flux limit 
Fth,k = Lth,k/^T^Dl, Lcut is the high luminosity cut-off of 
the XLF. As it is obvious from eQ. (|16|l and (|17p . the so com- 
puted /(L) is independent from SFRfc. 

To compute ^fe we used a two-step iterative procedure. 
At the first step w e assum ed 7 = 1.6, L cut = 2.1 • 10*" 
erg/s (|Grimm, Gilf anov, & SunyaevI |2003| ) and computed 
the XLF using eas.(IT6|l and (fT7|l . The Lth,fe was set to 1.5 
times the completeness limit for each galaxy. The obtained 
luminosity function was fit with a broken power law 



dL38 ^ 




71 r — f2 
^S8 ' 



Lag < Li, 



Lb ^ Las ^ Lc 



(18) 



where Lag = Lx/10 erg/s, L5 is the break luminosity, 
SFR is in units of A/0 yr~^ and ^ is the average normaliza- 
tion. The fit was done on unbinned data using maximum- 
likelihood (ML) method, the predicted contribution of re- 
solved CXB sources was included in the model as described 
below. Using this model and its best-fit parameters we cal- 
culated more accurate values of (,k and repeated the proce- 
dure. We obtained the following best-fit values for the XLF 
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Figure 5. The luminosity functions obtained by combining data of all the galaxies of the primary sample, normalized by SFR. Left: The 
top histogram with error bars (open circles, red) shows the luminosity distribution of all compact sources. The two histograms peaked 
near \og{Lx) ~ 37 are predicted contributions of resolved CXB (blue) and LMXB (ma genta, dashed) sources. The latter was computed 
using the average scaling relations for nearby early type galaxies from [Gilfanovl ||2004| ) . irrespective of their age. Apparently, it grossly 
overestimates the LMXB numbers, therefore no attempt to subtracts their contribution was made. The histogram marked with solid 
circles (black) shows the luminosity distribution of compact sources with the CXB contribution subtracted. It is our best approximation 
to the average HMXB XLF in nearby galaxies. The histogram marked "SFR" (green) shows the total SFR (in unit of lOOM0yr~^) of 
all galaxies contributing to a given luminosity bin in the XLF. Right: The CXB subtracted luminosity distribution from the left panel. 
The shaded area indicates the range of uncertainties corresponding to 40% uncertainty in the CXB normalization. The solid line on both 
panels shows a power law with parameters according to eq. lllSll . 



parameters: 71 = L58 ± 0.02, 72 = 2.73+0 54, = IIOI3I, 
^ = L49 ± 0.07. The value of the high-luminosity cut-off 
was fixed at Lcut = 5 x 10'^, the result being virtually 
independent on this choice. The final luminosity distribu- 
tion is shown in Fig. O The obtained XLF and its pa- 
rameters are remarkably simi l ar to the one determined by 
iGrimm. Gilfanov. fc SunvaevI (|2003l ). Correspondingly, the 
second step of this procedure does not has any noticeable 
effect. 

The statistical uncertainty of the slope 71 of the aver- 
age XLF, quoted above is rather small, thanks to the large 
number of sources involved in its determination. The rms of 
individual slopes around their average value is much bigger, 
and is investigated in the section VTM 

7.2 High-luminosity break 

The best fit XLF model has a break at Lx « 1.1 x 10*° erg/s. 
To access its significance we compare observed number of 
bright sources with the predictions of the single slope power 
law model. The best fit parameters of the latter are: 7 = 
1.60 ± 0.02, ^ = 1.49 ± 0.07. This model predicts 21.1 
sources above Lx = 10**° erg/s, whereas we detected 11 
sources. The discrepancy between the data and the model is 
~ 2.5ct, without account for the arbitrariness in the choice 
of the threshold luminosity, i.e. is marginally significant. 



The error contours for the break luminosity and the 
XLF slope after the break are shown in Fig.[71 demonstrating 
that the break parameters are not very well constrained, 
which is not surprising given its low statistical significance. 

7.3 Contribution of LMXBs and CXB sources 

Expected contributions of resolved CXB and LMXB sources 
to the XLF are shown in Fig. [S] They were computed sum- 
ming up individual predictions for all galaxies contributing 
sources to the given luminosity bin and then renormalizing 
the total in the same way as the XLF itself. For LMXBs 
we used the average s c aling relation for nearby early type 
galaxies from lCilfanwl (120041} . for CXB sources we used the 
log A'' — log S determined bv lGeorgakakis et al.l l|2008l ') . 

Both predictions bear uncertainties which may affect 
the resulting HMXB XLF. In the case of the CXB contri- 
bution, the uncertainties are caused by angular fluctuations 
of the density of background AGN and are likely to be re- 
duced as a result of averaging over a rather large number of 
galaxies involved in the calculation of XLF. However, their 
final amplitude is difficult to estimate, we therefore illustrate 
their possible effect by showing the range of XLF variations 
corresponding to 40% variations in the CXB normalization 
(shaded area in the right-hand panel in Fig[5]). As clear from 
the FigO the CXB contribution is significant and its uncer- 
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Figure 6. The CXB subtracted XLF of compact sources, multi- 
plied by LI 'S. The shaded are has same meaning as in Fig[5] 
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Figure 7. The error contours for the XLF break parameters: 
break luminosity L5 and the power law index after the break 72. 
The cross indicates the best fit values. 



tainties are relatively important in the log(L) ~ 36 — 38.5 lu- 
minosity range. It is also clear from the figure that the CXB 
contribution is predicted with a reasonable accuracy - the 
amplitude of the feature in the total XLF at log(L) ~ 36 — 37 
is consistent with the expected CXB contribution. 

Estimation of the LMXB contribution is less straight- 
forward. There is an evidence that the LMXB population 
may have a rather strong age dependence, being smaller for 



younger ga laxies jBogdan fc Gilfanov! [20T0I ). As the scaling 
relation of lCilfanovl ( 2004 ) was computed for a sample of 
nearby early type galaxies irrespective of their age, it may 
be not directly applicable to galaxies with significant ongo- 
ing star-formation. Indeed, the calculation based on results 



of ICilfanovl (|2004f ) apparently over-predicts LMXB contri- 
bution, probably by the factor of a few - its subtraction 
from the XLF would result in negative values (Fig[5}. For 
this reason we decided to subtract only the contribution 
of CXB sources, until detailed information on scaling re- 
lations for LMXBs in late type galaxies becomes available. 
The part of the XLF affected by this uncertainty is in the 
log(Lx) ~ 36.5 — 38.5 luminosity range. 

We note that the bright end of the XLF is virtually free 
of both contaminating components - the vast majority of 



sources brighter than 
mass X-ray binaries. 



10 erg/s in our sample are high 



7.4 XLFs of individual galaxies 

We fit XLFs of individual galaxies with a single slope power 
law with a high luminosity cut-off 



dN _ 



L < I/cut = 10 ^ erg/s 



(19) 



The cut-off luminosity was fixed at the value, exceed- 
ing the luminosity of the brightest source in the sample. 
Only galaxies having more than 5 sources above the detec- 
tion threshold were analyzed, the detection threshold be- 
ing defined as 1.5 times the 80% completeness limit of a 
gi ven galaxy. The CXB con tribution, computed according 
to ICeorgakakis et al.l l|2008l ). was included as a fixed com- 
ponent of the model. The fit was done with the unbinned 
data, using the Maximum Likelihood method. 

The best-fit values of slope are shown in Fig|8] The 
average (computed using method) is {7) = 1.59 with 
rms = 0.25. The = 16.8 for 24 de grees of freedom in- 
dicates that the data is consistent with the 7 being same 
in all galaxies. The two galaxies deviating from the mean 
by ~ 2a are NGC 3079 (7 = 2.62j:E;-^;^) and NGC 3310 
(7 = 1.30l°i^). 

There is some evidence that the XLF slope depends 
on the inclination of the galaxy (right panel in Fig[8]), but 
it is not statistically significant. The x^ = 14-2 for a lin- 
ear model 7 = 70 + fc X cos(i) with the best fit parameters 
of 70 = 1.71 ± 0.09 and k = -0.22 ± 0.14. According to 
the F-ratio test such a reduction of x^ corresponds to a 
« 2(7 fluctuation. We note however, that the character of 
the correlation (steeper XLF slopes for more edge-on galax- 
ies) is in accord with the theoretical expectation. Indeed, as 
the brighter sources (black holes in the soft spectral state) 
tend to have softer spectra, they are more subject to absorp- 
tion which is stronger in more edge-on galaxies. This makes 
the bright end of the XLF to be suppressed in more edge- 
on galaxies, i.e. their XLFs are expected to be steeper, as 
observed. No other physically meaningful correlations were 
seen in the data. 

Dependence of the XLF normalization A on the SFR is 
shown in Fig[9l The x^ fit to the logarithm of the normaliza- 
tion in the form log^ = a -I- /3 log SFR yields /3 = 0.73 ±0.05 
with very large value of the x^ = 99.2 for 23 d.o.f and the 
rms = 0.32 dex. Forcing a linear relation A — £^ x SFR we 
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Figure 8. Left: XLF slopes for individual galaxies plotted against the SFR (left) and inclination (right) of the galaxy. The horizontal 
solid line shows the mean value of (7) ~ 1.6. 



obtained ^ = 1.88 with significantly larger = 130.4 for 
24 d.o.f., but with rms = 0.34 dex, just marginally worse 
than obtained in the non-linear fit. Both best-fit models are 
shown in Fig[9] 

The relation between A and SFR can be expressed in 
terms of the relation between the number of sources above 
a certain luminosity limit and the star-formation rate: 

iVxRB(L > lO^^ergs"') = 3.22 x SFR(M0yr~') (20) 

For our sample of galaxies this relation predicts the number 
of sources with the accuracy of rms = 0.34 dex i.e. within a 
factor of ~ 2 (rms). 

The /? 7^ 1 obtained above implies a possible non-linear 
slope of the relation between XLF normalization (i.e. the 
number of sources) and SFR. However, neither linear nor 
non-linear model describes the data adequately {y^ /d.o.f . ~ 
4 — 5), suggesting that additional parameter(s) play a signifi- 
cant role in this relation. Therefore it is not clear how much 
weight should be given to the non-linear best fit slope of 
the A'x — SFR relation obtained above. It may, for example, 
be a result of application of an inadequate or incomplete 
model to a particular selection of galaxies. Alternatively, 
it may point at a genuine non-linear relation between the 
number of sources and SFR. A clue may be provided by 
accurate characterization of the factors which can poten- 
tially affect the number of HMXBs in star-forming galax- 
ies, such as metallicity and/or recent star- formation history. 
With the presently available data, it does not seems to be 
possible t o discriminate between these two pos sibilities. We 
note that iGrimm. Gilfanov. fc SunvaevI (|2003l ). based on a 
much smaller sample of galaxies, obtained a linear A^x — SFR 
relation. Finally, in interpreting the Ax — SFR and Lx — SFR 
relations one should keep in mind that for the XLF slope of 




SFR [M^yr-i] 

Figure 9. The XLF normalization versus SFR. The solid line 
shows the best-fit linear relation, dashed line - the best fit relation 
in the form A = g X SFR*^. See text for details. 



~ 1.6 the number of sources is defined by faint sources, near 
the detection limit whereas their total luminosity - by the 
(sometimes just a few) brightest sources. 
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8 COLLECTIVE LUMINOSITY OF HMXBS 

The coUective X-ray luminosities of X-ray binaries in galax- 
ies from our primary sample were computed as the sum of lu- 
minosities of point sources detected above the corresponding 
sensitivity limit Liim- The latter was defined as the luminos- 
ity level corresponding to a value of the incompleteness func- 
tion K{L) = 0.6, in the regions selected as described in Sect. 
|4] As Chandra observations of different galaxies achieved dif- 
ferent sensitivity, these luminosities need to be transformed 
to the same threshold, which was chosen at Lx ~ 10^® erg/s, 
i.e. the luminosities presented below are an estimate of the 
combined luminosity of X-ray binaries brighter that this 
threshold. The luminosity of "missing" unresolved sources 
was computed by integrating the luminosity function which 
differential slope was fixed at the average value of 7 = 1.58, 
determined in section [T] The normalization was determined 
for each galaxy individually, from the number of sources de- 
tected above its sensitivity limit, corrected for incomplete- 
ness and subtracting the CXB contribution as described in 
section[71 Finally the luminosities were corrected for the con- 
tribution of the background AGN to the resolved sources: 



fXRB 

^0.5-8 



Jl036 



E 



K{Li) 



AnD'FcKB 



(21) 

where -Fcxb is the predicted total flux of CXB sources 
brighter than the fiux limit correspo nding to I/Um computed 
from the log A*' — log S of .Georgakaki s et al.l (|2008l ) and D is 
the distance to the galaxy. The CXB contribution did not 
exceed ~ 35 — 40% of the total luminosity of galaxies, in ac- 
cord with the procedure used to define the analysis regions 
(Section 14. 3|) . For the reasons described in Section [7.31 we 
did not attempt to subtract the contribution of LMXBs. 

The obtained -Lxrb — SFR relation is shown in the left- 
hand panel of Fig. 1101 We also show in this plot the data for 
unresolved galaxies from the high-SFR sample. The latter, 
however were not used for the fits because of the contribution 
of unresolved emission, as discussed in section A least- 
squares fit to the data with the relation log Lx = log K + 
P log SFR yielded a value of slope consistent with unity: /3 = 
1.01 ± 0.11. We therefore fixed the slope at the unity and 
repeated the fit to obtain the best-fitting linear relation: 

iaf-8kcv(ergs-^) = 2.61 x 10^^ SFR (Mq yr'^) (22) 

The dispersion around this relation is ct = 0.43 dex. 

To facilitate comparison with previous results, we also 
studied the Lx — Lm, relation (right-hand panel of Fig. llOp . 
The linear fit of this relation gives the following result: 

ia5-8kev(ergs"^) = 1.75 x 10"" Lir (ergs"^) (23) 

where Lir is the total (8—1000 \im) infrared luminosity. The 
dispersion around the best-fit is u = 0.54 dex, i.e. notably 
worse than for the Lx — SFR relation with SFR based on 
combined UV-I-IR data. 



8.1 Unresolved galaxies 

As it is obvious from Fig llOl unresolved galaxies are dis- 
tributed along the extension of the Lx — SFR relation for 
resolved galaxies, albeit somewhat higher. Quantitatively, a 
linear fit to the Lx — SFR relation using unresolved galaxies 



only yielded the scale factor of ~ 3.7 x 10^^ erg/s/(M0 /yr), 
i.e. by a factor of ~ 1.5 higher than in eq. (|22p . This re- 
sult is to be expected because the luminosities of unre- 
solved galaxies, by definition, include the unresolved emis- 
sion component, which is excluded in the case of the galaxies 
from our main sample. Although the detailed study of un- 
resolved emission is beyond the scope of this paper and will 
be fully addressed in the Paper II (in preparation), we made 
a preliminary estimate of its amplitude in several galaxies 
from the primary sample, uniformly distributed along the 
SFR axis. We found that unresolved emission contributes 
~ 20 — 60% to the total X-ray luminosity of these galax- 
ies. This is consistent with the higher normalization in the 
scaling relation for unresolved galaxies, obtained above. 

8.2 Comparison with previous studies 

To do the comparison, we converted the results of the earlier 
studies to the 0.5-8 keV band assuming an absorbed power 
law with r = 2.0 and = 3 x 10^^ cm"^ 

The first systematic studies of the Lx— SFR scaling rela- 
tions based on the data o f mod er n X-ray satellites were per- 
forme d by iBauer et al.l ||2002|). Ranalli. Comastri. fc Settil 
( l2003h and lGrimm. Gilfanov. fc SunvaevI {200^ . 

In, pro bably the first, Ch andra-based publication on 
this subject, iBauer et al.l (I2OO2I ) studied the correlation be- 
tween X-ray and 1.4 GHz luminosity of faint X-ray sources 
in the Hubble Deep Filed North. Using their eq.(5) we find 
that in the SFR range of interest their Lx/SFR ratio exceeds 
by a factor of ~ 3 — 4 the value obtained in the present work. 
As this is also inconsistent by about same factor with the 
results of many other studies, we will not attempt to inves- 
tigate the origin of this discrepancy. We note however, that 
the scale factor of their 1.4 GHz -based SFR determination 
is similar to the one used in this paper, which produced 
con sistent results for the CDF-N gala xies . 

iRanaUi. Comastri. fc Settil l|2003l ) obtained tight corre- 
lation between X-ray, radio and FRI luminosities for a sam- 
ple of nearby galaxies using X-ray data from ASCA and 
BeppoSAX satellites. Using their eqs (8) and (12), convert- 
ing X-ray luminosities from 0.5-10.0 to 0.5-8.0 keV band (a 
factor of 1.11) and FIR (42.5— 122.5 vim) luminosity to total 
IR (8—1000 vim) band using a conversion factor of 1.7, we ob- 
tain the relation Lo.5-8keV ~ 2.4- 10"'* x Lir. The scale fac- 
tor in this relation is by « 1.36 l arger than in eq.(|23p. which 
is con sistent with the fact that iRanalli. Comastri. fc Settil 
(|2003h analyzed the total luminosities of galaxies, includ- 
ing unresolved emission. We conclude that their Lx — Lir 
relation is consistent with the one derived in this paper. 

To derive Lx — SFR relation we use their eqs. (14) and 
(15) and obtain Lo.s-skcV ~ 8.6-10^^ x SFR(> 5M0) where 
SFR(> 5M0) is the formation rate of massive stars. For 
Salpeter IMF SFR(0.1 - IOOMq) = 5.5 ■ SFR(> 5M0), we 
thus obtain Lo.5-8koV ~ 1-6 ■ lO^'' x SFR(0.1 - IOOM0). 
This is close but not identical to the relation derived in the 
present paper, ea. ((22)) . especially considering that it applies 
to the total X-ray luminosity of the galaxy, not only bright 
X-ray binaries. As the underlying scaling relations between 
X-ray and IR luminosities are fully consistent with our re- 
sults, the factor of ~ 2 — 2.5 discrepancy in the Lx — SFR 
relation must be due the different method of the SFR esti- 
mation used in the two studies. 
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Figure 10. The Lx ~ SFR {left) and Lx ~ ^IR ('"'ff'**) relations. In both panels the squares are galaxies from the resolved 
sample, the triangles and filled circles are galaxies from the unresolved high-SFR sample. The solid lines show the best-fitting 
scaling relations, obtained using only resolved galaxies and given by eqs. i'l'll and I I23I I respectively. 



used Chandra ob- 
servations of nearby (i.e. spatially resolved) star-forming 
galaxies and derived a relation between luminosity of com- 
pact sources and SFR. Using their eq.(21) and converting 
the luminosity to 0.5-8 keV band (conversion factor of 1.28) 
we obtain: Lo.5-8kov « 8 .5-10^^ xSFR(M > 5Ma),i .e. iden- 
tical to that obtained by 
IShtvkovskiv fc Gilfanov! ( 



Ranalli. Comastri. fc Settl (|2003h . 
20051 ) recomputed the scale in this 



relation with account for the updated calibration of IR and 
radio SFR indicators (Bell 2003). Using their scale factor we 
derive Lo.5-8kev ~ 2.8 ■ 10^^ x SFR(0.1 - IOOMq), in good 
agreement with eq. (|22|) . 

As a historical note we mention that the 
SFR values de termined from the IR-b as ed c alibra- 
tion used by Ranalli. Comastri^ fc Settil (|2003h and 



iGrimm. Gilfanov. fc SunvaevI ( 20031 ) were incorrectly at- 
tributed by the authors to the formation rate of massive 
stars. This was a result of some confusion that existed 
in the SFR determination-related literature at the time 
of those publications and, in particular, by inconsistent 
definitions of radio- and IR-based SFR proxies. This error 
could potentially result in a factor of ~ 5.5 overestimate 
of the scale in the Lx — SFR relation. However, it was 
compensated by the usage of FI R instead of total IR 
luminosity in the lKennicufr3 l| 19981 ) relation and the scale 
factors somewhat different from the presently accepted 
values. As a result, the final Lx — SFR relations derived 
by these authors are in a reasonable agreement with the 
scaling relations obtained in this paper. 

Among more recent public ations on this subj ect, w e 
compare with the resuhs of IPersic fc~ R cphacli J2003) , 
iKaaret fc Alonso-Herrerd (|2008l ) and lLehmer et all (|201G[ ). 

IPersic fc Rephaelil ^M) 

investigated the relation be- 
tween the collective hard X-ray luminosity of young point 
sources and the SFR derived from infrared luminosity. Using 



their eq.(9) and converting the 2-10 keV luminosity to the 
0.5-8 keV band (a factor of 1.28) we obtain: Lo.s-skcV ~ 
0.96 • 10^^ X SFR(0.1 - IOOM0). This scale factor is signifi- 
cantly, by a factor of ~ 2 — 2.5 lower than derived in all other 
studies, including this one. The reason for this discrepancy 
is, in our opinion, in an overestimated contribution of the old 
pop ulation to the X-ray emiss i on of s tar-forming galaxies. 

iKaaret fc Alonso-Herrerol (j2008l ) examined the X-ray 
point-source population of two galaxies with high SFRs 
(NGC 4194 and NGC 7541) for which they found Lx/Lm 
ratio consistent with our average relation, eq. (|23|l . However, 
their Lx/SFR ratio is by a factor of ~ 2 larger than ours, 
due to a different SFR calibration used. 

Based on Chandra observations of nearby luminous in- 
frared galaxies (LIRGs) Lehmcr ct al. (2010) studied the 
relation between 2-10 keV luminosity and SFR. To com- 
pare results, we convert their 2 — 8 keV luminosity to 
the 0.5 — 8 keV band and obtain the scale factor for the 
linear Lx — SFR relation (first line in their Table 4) of 
Lx/SFR ^ 2.2 • 10^^ergs"^(Mgyr-^)-\ that is in a good 
agreement with our result. 



8.3 Effects of statistics of small numbers in the 
Lx - SFR relation 

For the power law luminosity function with a slope of 1.6, 
the Lx — SFR relation is expected to be non-linear in its low- 
SFR part, due to the effects of statistics of small numbers 
(Gilfanov, Grimm, fc Sunyaev 2004a, 2004b). These effects 
are important when considering the collective luminosity of 
a population of discrete sources and are a consequence of an 
asymmetric probability distribution for the sum of luminosi- 
ties of compact sources in the limit of their small number. 
In the case of HMXBs in star-forming galaxies, they make 
the Lx — SFR relation to appear non-linear in the low SFR 
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limit. This effect is purely statistical in its nature, in the 
sense that it affects the most likely value of the luminosity 
of a randomly chosen galaxy, whereas the relation between 
average luminosity of many galaxies and SFR is linear in the 
entire range of star-formation rates. This is achieved due to 
extended tail of the probability distribution j3(Ltot) towards 
large Ltot- In the sample of galaxies from Gilfanov, Grimm, 
& Sunyaev (2004a, 2004b) such a non-linear dependence was 
clearly observed at SFR < 5 Mq /yr. 

Despite nearly identical XLF of compact sources ob- 
tained in this paper, the Lx — SFR relation appears to be lin- 
ear in the entire SFR regime. In an attempt to find the rea- 
son of such behavior, we investigated the distribution of the 
luminosity of the brightest source in a galaxy and found that 
for a fraction of low-SFR galaxies it is higher than predicted 
based on the average XLF ( the predicted value computed 
accord ing to formulae from iGilfanov. Grimm, fc SunvaevI 
l|2004f) ). As XLF slopes of low-SFR galaxies are consistent 
with the rest of the sample (Fig|8]) and the contribution 
of background AGN is negligible in the bright end of the 
XLF, we conclude that this may be a consequence of the 
observer bias, introduced by selection of galaxies - targets 
for Chandra observations (see Sect. [2^ . For example, the 
proposers may have predominantly selected galaxies which 
were known to have enhanced level of X-ray emission, e.g. 
based on the ROSAT all-sky survey. This issue will be in- 
vestigated in more detail in a follow-up publication. 

8.4 Dispersion in the Lxrb — SFR relation 

Despite of the special effort to produce a sample with a 
homogeneous set of multi-wavelength measurements and to 
minimize contamination by CXB sources and LMXBs, the 
resulting -Lxrb — SFR relation bears a rather large disper- 
sion of ~ 0.4 dex rms. In order to investigate the origin of 
this dispersion we searched for correlation of the Lxrb /SFR 
ratio with different quantities characterizing importance of 
various potential contaminating factors. These included the 
sensitivity limit of Chandra data Lum, contribution of CXB 
and LMXB sources characterized by the ratio of the solid 
angle of the studied region and enclosed stellar mass to the 
SFR, dust attenuation characterized by the Lir/Luv ratio 
and inclination of the galaxy i. The results are shown in Fig. 

El 

We analyzed these data using the Spearman's rank cor- 
relation test and found no statistically significant correla- 
tions. The results of these tests are presented below with 
rs being the correlation coefficient and P being the prob- 
ability under the null hypothesis that the variables are un- 
related. We found no correlation with the point source de- 
tection sensitivity {rs ~ 0.33, P = 8%). Therefore we can 
exclude that the extrapolation of X-ray luminosity above 
10'^^ erg/s is the main reason for the dispersion around 
the Lxrb — SFR relation. There is no evidence for strong 
contamination by CXB sources (rs — —0.03, P — 86%) 
or LMXBs (rs = -0.03, P = 88%). This proofs that the 
spatial analysis, performed in order to minimize the contri- 
bution of CXB and LMXB sources to the compact source 
population, yielded sensible results. It is further confirmed 
by repeating the previous test on the data without CXB 
subtraction: rs = —0.01, P = 96%. We used the ratio of 
IR to NUV luminosity Lir/Lnuv to characterize the dust 
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SFR,„, [Mg yr-'] 

Figure 12. Comparison of the SFR estimations obtained by dif- 
ferent methods. The x-axis is the SFR used in the paper, com- 
puted as described in Section |6] The y-axis values are computed 
from the UV l uminosity c orrected for the dust attenuation using 
the method of lBuat et al.l 12005 ;^ . The straight line is y = x rela- 
tions. The dispersion of points around this relation is rms = 0.12. 



attenuation (higher the ratio, higher the absorption) and 
found virtually no correlation (rs ~ — 0.35, P = 7%). Nei- 
ther Lx/SFR correlates with the inclinations of galaxies 
(rs = 0.09, P = 65%). To verify this approach, we ana- 
lyzed the dependence between Lir/Lnuv ratio and incli- 
nation of the galaxy and found a strong anti-correlation 
(rs ~ —0.59, P = 0.2%), as expected. We note that in order 
to suppress the 0.5 — 8 keV luminosity by a factor of 10 a 
large column density of ~ 5 x 10^'' cm~^ is needed. 

We also tried to probe the amplitude of intrinsic uncer- 
tainties of the SFR determination. To this end we computed 
the UV based SFR estimates corrected fo r the a ttenuation 
effects following the method o f iBuat et al.l(|2005l ). and com- 
pared with the SFR values used in this paper (Fig |12|) . We 
found a good consistency between the SFR estimations ob- 
tained by the two methods with rms — 0.12 dex, i.e. three 
times smaller than the rms of the Lx — SFR relation. This 
suggests that uncertainties of the SFR determination are not 
the main source of the scatter in the Lx — SFR relation. 



9 DISCUSSION 

We have found rather large scatter around the average 
Lx — SFR relation rms ^ 0.4 dex (FigllOp. As discussed 
in the Section 18.41 the scatter is unlikely to be caused by 
any of the obvious contaminating factors, such as difference 
in the sensitivity limits of observations of different galaxies, 
pollution by CXB sources or LMXBs. This proofs that the 
criteria used for selecting the sample and the method for 
choosing the regions for the analysis yielded the desired re- 
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Figure 11. Dependence of the log(Lx/SFR,) on quantities cliaracterizing different potentially contaminating factors: the point source 
detection sensitivity of Chandra observations (a), the fractional contribution of CXB sources (b) and LMXBs (c), the dust attenuation 
characterized by Lir/Lnuv ratio (d) and the inclination of the galaxy. The last panel (e) shows the relation between the dust attenuation 
and inclination. Each point represents a galaxy from the primary sample. The crosses show average values in bins; the bin extent is 
indicated by the horizontal error bar, the vertical error bar shows the rms of points with respect to the average value. The contribution 
of CXB sources is subtracted. See Sect 18.41 for details. 



suit. On the other hand it suggests that the observed scatter 
has physical origin. 

Although similar scatter is found in the SFR- 
dependence of the XLF normalization, rms ~ 0.34 dex 
(Fig[9]), there is no precise one-to-one correspondence be- 
tween the deviations of these two quantities from their best 
fit relations, as illustrated by Fig llSI At the same time the 
XLF slopes appear to be consistent, within statistical er- 
rors, with all having the same (or close) value (FiglSl [T4)) . 
These two facts do not contradict to each other. Indeed, 
for a power law with 7 = 1.6 the XLF normalization and, 
largely, its slope are determined by the more numerous faint 
sources, whereas the total luminosity is determined by the 
(often few) brightest sources in the galaxy. Thus, we con- 
clude that the scatter in the Lx — SFR relation is caused by 
both variations in the XLF normalization and fluctuations 
in the number and luminosity of brightest sources in the 
galaxy (cf. Fig |14|l , the variations of the XLF slopes playing 
less significant role. 

Metallicity is one of the primary quantities that may af- 
fect the population of HMXBs, as it exerts an e qually great 
effect on HMXB creation and emission (e.g. IClark et al.l 
1 19781 ). It has a rather substantial impact on the strength of 



stellar winds, thus on the luminosity of individual HMXBs, 
making it larger for larger amount of metals present in the 
companion star. Star formation at lower metallicity is ex- 
pected to produce more massive stars as a result of the 
low cooling rate. This effect would result in a higher num- 
ber of compact objects in binary systems. This simple the- 
oretical argument is supported by the recently found evi- 
dence that the number of ultra-luminous X-ray sources is 
inversely proportiona l to the metallicity of their host galaxy 
(|Zampieri fc Roberts! l2009l : iMapeUi et al] I2OI0I ). However, 
this is in contrasts w ith population synthesis calculations of 
iLinden et all (|2010l ). who found that the existence of mas- 
sive single star BHs at low metallicities does not directly 
correlate to the existence of more massive compact black 
holes in binary systems. Despite this, they find the same 
inverse relationship between the number of ultra- luminous 
X-Ray sources and metallicity, described in iMapelli et al.l 
(2010). In an attempt to advance this issue observationally, 
we collected from the literature a rather heterogeneous set 
of the metallicity measurements for a relatively small subset 
of galaxies, but could not find any meaningful trends in the 
data. This, however, should not be considered as an ulti- 
mate result, more systematic studies are obviously needed 
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Figure 13. The XLF normalization vs. total luminosity of X- 
ray binaries. Both quantities are normalized to the SFR with 
the normalization coefficients equal to the best fit values in the 
corresponding scaling relations. Only galaxies having 5 sources 
above their threshold luminosities are shown. This plot illustrates 
the lack of one-to-one correspondence between the deviations 
from the average dependences in the number of sources and in 
the total luminosity. 



in order to make a reliable conclusion regarding the role of 
metallicity (work in progress). 

The synthetic XLF produced by iLinden et all (|2010l ') 
has a high luminosity turnover at approximately 10*° erg/s, 
close to the observed value, though the slope of the latter 
XLF is somewhat harder (1.16) than observed . The p ossible 
reason for the latter may be that 'Linden et al. ' (2010) do not 
consider variable HMXBs such as Be-HMXBs which provide 
much of the lower-luminosity population (T. Linden, private 
communication) . 

Another important factor is the recent star-formation 
history of the host galaxy. Indeed, as more massive stars 
have shorter life times, and, on the other hand, have stronger 
winds one may expect that the very luminous HMXBs are 
short-lived. This effect may result in the dependence of the 
X-ray luminosity of the galaxy on the recent star-formation 
activity, which is not characterized precisely the SFR value 
determined from the IR and/or UV data, thus introducing 
scatter in the Lx — SFR relation. 

In addition to the more fundamental reasons discussed 
above, variability of the brightest X-ray sources may also 
contribute to the scatter in the Lx — SFR relation. The 
X-ray sources with logLx ^ 39 (usually classified as ultra- 
luminous X-ray sources) likely harbor a black hole and may 
be subject to the disk instability resulting in transient events 
characteriz ed by luminosity variation s of several orders of 
magnitude (|Tanaka fc ShibazakI 19961 ). In addition, spectral 
state changes may lead to the luminosity changes from the 
factor of a few upto ~order of magnitude. Their effect on 
the total luminosity will be especially important in galaxies 



with SF R < 10 M(7)/yr, containing only a few such bright 
sources (|Gilfanov. Grimm, fc SunvaevlbOO^ ). 

9.1 XLF of high-mass X-ray binaries and ULXs. 

Based on a sample of ~ 700 compact sources, we have 
produced the average luminosity distribution of HMXBs in 
galaxies (Fig[5])- This is a significant improvement in the 
st atistical accuracy as compared with the previous result 
of [g rimm. Gilfanov. fc SunvaevI (|2003h . We have also im- 
proved in terms of the errors introduced by the distance 
uncertainties, thanks to the systematic compilation of dis- 
tances for galaxies which became available in the NED. 
Thirdly, we made a special effort to minimize and accu- 
rately subtract the contribution of background AGN. Al- 
though we could not in the same way remove the contribu- 
tion of LMXBs, because of the strong age dependence of the 
LMXB-stellar mass scaling relations, the method of con- 
struction makes us confident that our sample is relatively 
free from the LMXB contamination. This is definitely true 
in the high luminosity end, logLx > 38.5 — 39. 

Despite all these improvements our average HMXB 
XLF is entirely consistent with the one obtained by 
iGrimm. Gilfanov. fc SunvaevI l|2003l ). Indeed, the values of 
slope agree to the accuracy of ~ 0.03 and the values of 
the high luminosity break are consistent within (admittedly 
rather large) s tatistical uncertainties. 

Similar to [Grimm. Gilfanov. fc SunvaevI (j2003l ). we did 
not find any statistically significant features in the XLF near 
the critical Eddington luminosity of a neutron star or a stel- 
lar mass black hole. The XLF follows a single slope power 
in the entire luminosity range, upto logLx ~ 40, where it 
breaks. Taken at the face value, the break corresponds to 
the Eddington luminosity of a ~ 50 — 100 Mq black hole. 
We note that the statistical accuracy of the XLF, grouped 
into broad, AlogL ~ 0.5, luminosity bins, is < 13% and 
< 22% in the log Lx ~ 38 - 39 and log Lx ~ 39 - 40 lumi- 
nosity r anges (Fig|6]). This further strength ens the propo- 
sition of[G rimm. Gilfanov. fc SunvaevI (|2003l ) that the bulk 
of the population of what is called the ultra-luminous X-ray 
sources (conventionally defined as sources with log Lx > 39) 
are a high luminosity end of the population of the ordi- 
nary HMXBs and most likely harbor stellar mass black holes 
rather than exotic intermediate mass objects. The very high 
luminosity end of the distribution, log Lx > 40, may, how- 
ever, be associated with the intermediate mass black holes. 

The absence of features in the luminosity function of 
HMXBs is still puzzling. The inaccuracies in distances to 
the galaxies in our sample may still contribute to some 
smearing of the break at the NS Eddington limit, but 
th eir role must be significantly smalle r than in the XLF 
of [g nmm. Gilfanov. fc SunvaevI l|2003h . Indeed, the aver- 
age distance uncertainty for the galaxies from our sample is 
~ 15% (Table 111 which translates to luminosity uncertainty 
of A log(Lx) ~ 0.12. This corresponds to the half of the bin 
width for the XLF shown in Fig|5] (and < 1/4 of that for the 
one shown in Fig[5]) and can not have any major effect on the 
XLF shape. Taken at the face value, this result may suggest 
that there is no sharp limit on the luminosity of an accreting 
compact o bject at the Eddington value, as discu ssed, for ex- 
ample, bv iGrimm. Gilfanov. fc SunvaevI (|200l [20031 This 
does not explain, however, why the nearly an order of mag- 
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Figure 14. Several representative XLFs, normalized to the SFR. 
In the order of increasing luminosity of the brightest sources: 
MlOl, NGC 3079, M51, Antennae, NGC 3310. The solid line 
shows a power law with average normalization of A = 1.88, 
slope 7 = 1.6 and a cut-off at Lcut = 10*" erg/s. The XLFs 
of NGC 3079 and NGC 3310 are, correspondingly, the steepest 
and the flattest in the sample. 



nitude difference in the specific frequencies of neutron stars 
and black holes does not reveal itself in the XLF. Indeed, 
in the entire luminosity range the XLF slope is consistent 
with the slope of the distribution of the mass transfer rates 
in HMXBs, p redicted from the standard Salpeter IMF for 
massive stars (|Postnovll2o"o3 ). 

Answers to these questions could be obtained if one 
could disentangle the contributions of NS and BH binaries 
and separate their XLFs. Direct measurements of the mass 
of the compact object in a large number of binaries in exter- 
nal galaxies is not feasible. Therefore indirect methods will 
have to be utilized, for example based on the X-ray spectral 
properties of the binaries. For example, one of such meth- 
ods, based on the X-ray colors, is being developed by Farrell 
et al (in preparation). 



10 IMPLICATIONS FOR THE THEORY OF 
BINARY EVOLUTION 

10.1 Specific frequency of X-ray bright compact 
objects in HMXBs 

With the knowledge of the relation between the number of 
HMXBs and SFR, we can estimate the fraction of compact 
objects that went through an X-ray active phase powered by 
accretion of matter from a massive donor star in a binary 
system. 

According to eq. psp . the number of HMXBs with lu- 
minosity higher than 10"^^ erg/s is: 

iVHMXB(> lO^^ergs"') 135 X SFR (24) 



On the other hand, the number of HMXBs at any given time 
is: 

A^HMXE ~ Nco ^ /x,fe T-x,fe ~ Nco /x tx (25) 

k 

where Nco ~ Nt{M > 8 Mq) is the birth rate of compact 
objects (both NS and BH), approximately equal to the birth 
rate of massive stars N^,{M > 8Mq). For a Salpeter IMF, 
we obtain: 

N,{M >8Mq) ^ 7.4 • 10"^ X SFR (26) 

The summation in eq. (|25p is done over different types of 
HMXBs (wind-fed and Roche lobe overflow systems with a 
supergiant companion, Be/X binaries etc), fx,k is the frac- 
tion of compact objects (computed with respect to their to- 
tal number in the galaxy) which become X-ray sources due 
to accretion in HMXBs of the fc— type, rx,fc is the duration 
of their X-ray active phase. The fx = X^fc /x,fc is the to- 
tal fraction of X-ray bright compact objects and fx is the 
average life time of X-ray sources: 

^ T fx A ^ ^ 

/ ./ X,A; 

From binary evolution calculations, rx ~ 10* yr for super- 
giant systems and rx ~ 10^ yr for Be/X binaries. Majority 
of known high-mass X-ray binaries are Be/X systems, espe- 
cially in Magellanic Clouds, where only a few supergiant 
syste ms are known ([Liu, van Paradiis. fc van den Heuvell 
l2006h . The fraction of the latter is larger in the Milky Way, 
where about 29 out of 114 HMXBs (are suspected to) have 
supergiant donors. If these proportions are representative for 
star- forming galaxies in general, the summation in eq.((25| 
and H27|) is dominated by Be/X systems and appropriate 
scaling for fx should be ~ 10"" yrs. 

Combining equations 1241 1251 and 1261 we obtain: 

This quantity represents the fraction of compact objects that 
on the time scale of ~ 10^ Myrs after their formation (e.g. 
IShtvkovskiv fc Gilfanovl 120071 ') become X-ray sources with 
Lx > 10^^ erg/s, powered by accretion from a massive donor 
star in a high-mass X-ray binary. 

It is interesting to note, that if the HMXB XLF contin- 
ues towards lower luminosities with ~ same or similar slope, 
it may be possible to constrain fx based on the specific fre- 
quency of X-ray sources in star-forming galaxies. Indeed, if 
XLF from ea. (|18[) can be extrapolated down to Lx ^ 10^* 
erg/s, as suggested by IShtvkovskiv fc Gilfanovl (|2005l ). the 
coefficient in eq. (|28p becomes 0.72. As fx can not exceed 
the binary fraction ftin ~ 0.5, the X-ray life-time of Be/X 
systems is fx < 0.07 (/bm/0.5) Myrs. 

10.2 Bright sources — Black hole systems 

The derivation of eq. (|28p assumed that fx and fx do not 
depend on the mass of the compact object which assump- 
tion is likely to be incorrect. However, as we used a rather 
low luminosity limit of 10^^ erg/s in ea. (|24|l . the A''hmxb 
is dominated by faint Be/X binaries, which harbor neutron 
stars. Therefore the fx derived above characterizes neutron 
stars and is likely to be insensitive to the contribution of 
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black holes. As contributions of BH and NS systems to the 
A^'hmxb in ea. (|24p can not be separated, we can not estimate 
/x for BH binaries of all luminosities separately. However, 
this can be done for bright sources, with luminosity signif- 
icantly exceeding the Eddington limit for the neutron star. 
Using 10^^ erg/s as the luminosity limit, we obtain: 



A''hmxb(> 10 ■ ergs 



; 0.49 X SFR 



(29) 



As we consider luminous sources, the Nco in ea. (l25p is now 
the birth rate of black holes, which for Salpeter IMF is: 

-3 



iV.(M > 25 Mo) 



1.4- 10"'' X SFR 



(30) 



The resulting fraction of stellar-mass BH that became lumi- 
nous X-ray sources in HMXB is: 



/bright 
X 



^ 3.5 ■ 10-^ X ' ™« 



(31) 



V lO"* yr 

In the latter equation we took into account that bright 
HMXBs shining near the Eddington limit typically contain 
supergiant donors, therefore their characteristic life times 
are shorter than for Be/X systems and are ~ 10* yrs. 



Finally we note that in the estimates in sections 110.11 
and ll0.2l we intentionally used Salpeter mass function rather 
than a more realistic IMF. Indeed, the UV and IR SFR 
proxies used in this paper rely on the emission from mas- 
sive, M > 5Mq, stars. The star-formation rates derived by 
these proxies are then converted to the full 0.1 — 1OOM0 
range assummg Salpeter IMF ('e.g. lKennic"ufr3ll998l ). There- 
fore, for the sake of consistency, same IMF should be used in 
the estimates of the numbers of compact objects in eq. (|26p 
and (|31|l . In the remaining part of this section a more cor- 
rect approach is to use a realistic IMF, correctly predicting 
the frequencie s of low mass stars, e.g. the one suggested by 
l|Kroupall200l ). 

10.3 Specific frequency of X-ray bright compact 
objects in LMXBs 

For comparison, we obtain a similar estimate for the frac- 
tion of compact objects that become X-ray sources pow- 
ered by accre tion from a low- mass donor star in an LMXB. 
According to iGilfanovl (|2004l '). the number of LMXBs with 
luminosity higher than 10'^^ erg/s is: 

iVLMXB(> lO^'^ergs-') ^ 50 x M,/{10'-°Mq) (32) 

Assuming stationarity: 



A^LMXB ~ Nco fx 



(33) 



Strictly speaking, stationarity assumption does not apply to 
bright LMXBs. Indeed, the life time of a source with Lx ~ 
10^^ erg/s (M ~ 10~* Mq/jt) powered by accretion from a 
~ 1 Mq star can not exceed ~ 10* yrs. However, as in the 
case of HMXBs, the A^lmxb in eq. (|32p is dominated by faint 
sources, for which the stationarity assumption may apply. 
Assuming a Kroupa IMF (slope of 7 = 1.3 for M < O.5M0, 
- y = 2.3 for 5 < M < IMq and 7 = 2.7 for M > IMq 
l|Kroupal2002D '). the number of neutron stars and black holes 
in such a galaxy would be: 

Nco ~ iV. (M > 8 Mq) ^ 5 ■ 10"^ X M, (34) 

Combining eqs. (|32p -(l34 p we obtain: 

/x ~ lO-'^ (35) 



According to the derivation, this number takes into account 
only LMXBs in bright state and ignores compact objects in 
transient systems which are in the quiescent state. 

It is interesting to note that the /x for LMXBs is by 
~ 5 orders of magnitude smaller than for HMXBs. This is 
another manifestation of the fact that LMXBs are extremely 
rare objects and may be explained by the high probability of 
disruption of the binary system with a low mass companion 
in the course of the supernova explosion (e.g. Verbunt & van 
den Heuvel 1995). 



10.4 Constrains on the mass-ratio distribution in 
binaries 

As not every compact object is in a binary and not every 
compact object in a binary has a suitable combination of bi- 
nary parameters to become an HMXB, fx derived in section 
llU.ll obevs the condition: 



/x</Mn(m2>5MQ) 



(36) 



where fbin(m2 > 5M0) is the fraction of binaries with com- 
pact object, having a massive companion. This inequality 
can be used obtain interesting constrains on the mass dis- 
tribution of the secondary in massive binaries. 

Firstly, we can exclude the possibility that the mass 
of the secondary in a binary system obeys a Kroupa mass 
function. Indeed, for Kroupa MF, the fraction of stars more 
massive than 5Mq is /(m > 5M0) « 5.5 • 10~^ by number, 
i.e. is by a factor of ~ 30 to ~ 100 smaller than fx, in a 
strong contradiction to ea. (l36|l . Considering a more general 
power law mass distribution tl}{m) oc m~^, and assuming 
fx ~ 0.2 we obtain from eq. (|36p a constrain 7 < 0.3, i.e. 
much flatter than the high mass end of Kroupa or Salpeter 
mass function. 

If the mass transfer in the course of the binary evolution 
does not significantly change the mass of the secondary, the 
above result strongly excludes the possibility that the masses 
of stars in a massive binary are drawn independently from 
the IMF. Assuming further that number of binaries obeys 
the distribution in the form 



dN oc ip{mi)dmi f{q)dq 



(37) 



where mi is the mass of the primary, with the mass dis- 
tribution tp{mi) and f{q) is the distribution of the mass 
ratio q — m-ijmx, we write for the fraction of system with 
the secondary mass m-i > 5Mq among all binaries with the 
primary mass mi > SAf© (i.e. those producing a compact 
object): 

draxMrax) C , f{q)dq 

fb.n{m2 > 5Mq) ^ 100 / (38) 

Jo.i ip{m)dm 

For the Kroupa IMF and /(g) = 1 the above equation gives 
fbin{'ni2 > 5M0) « 0.6, i.e. consistent with fx ~ 0.2. 



11 CONCLUSIONS 

Based on a homogeneous set of X-ray, infrared and ul- 
traviolet observations from Chandra, Spitzer and GALEX 
archives, we studied the properties of populations of high- 
mass X-ray binaries in a sample of 29 nearby star-forming 
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galaxies and their relation with the SFR. We took a spe- 
cial care to minimize the contamination by LMXBs, back- 
ground AGN and to control the incompleteness of the Chan- 
dra source lists. 

Our main findings can be summarized as follows: 

i) We have detected 1057 compact X-ray sources, of 
which ~ 300 are expected to be background AGN. The ma- 
jority of remaining ~ 700 sources are young systems associ- 
ated with star-formation in the host galaxy. Based on their 
high X-ray luminosities and analogy with the X-ray popula- 
tions in the Milky Way and few other very nearby galaxies, 
we conclude that they are high-mass X-ray binaries, pow- 
ered by accretion of matter from a massive donor star in a 
binary system. 

ii) The shape of their luminosity distribution is simi- 
lar in different galaxies with the power low indexes having 
rms — 0.25 with respect to the average value of « 1.6. 
The XLF normalization, on the contrary, show significantly 
larger dispersion with the rms = 0.34 dex around the 
A oc SFR law. 

iii) Combining the data of all galaxies we produce the 
average XLF. Despite its by far better statistical accuracy 
the obtai ned XLF is entirely consistent with the one pro- 
duced by iGrimm. Gilfanov. fc SunvaevI (|2003l ) based on a 
much smaller sample of galaxies and much fewer sources. 
The HMXB XLF has a slope of 7 « 1.6 in the broad lu- 
minosity range logLx ~ 35 — 40 and shows a moderately 
significant evidence for the high luminosity break or cut-off 
at log Lx ~ 40. Similar to Grimm et al, we did not find any 
statistically significant features at the Eddington luminosity 
limits of neutron stars or stellar mass black holes. 

iv) HMXBs are a good tracer of the recent star forma- 
tion activity in the host galaxy. Their collective luminosity 
is proportional to the SFR: 

ia5-8kov(ergs'^) = 2.61 x 10^^ SFR (Mq yr'^) (39) 

The rms of points around this relation is 0.4 dex. The ob- 
served dispersion is unlikely to be caused by any of the ob- 
vious contaminating factors such as CXB or LMXB sources 
and is likely to have a physical origin. 

v) The collective luminosity of X-ray binaries is well 
correlated with the infrared luminosity of the host galaxy: 

ij5-8kov(ergs-^) = 1.75 X 10"* Lm (erg s^^) (40) 

where Lir is the total (8—1000 [im) infrared luminosity. The 
dispersion around the best-fitting relation is a = 0.5 dex. 

vi) The number of bright HMXBs scales with the SFR 
of the host galaxy: 

iVxH,B(> 10^**ergs-^) = 3.22 x SFR(M0yr"') (41) 

with the dispersion of rms = 0.34 dex. 

vii) The fraction of compact objects which were X-ray 
bright at least once in their lifetime, powered by accretion 
from a massive star in a binary is fx ~ 0.2. This constrains 
the mass distribution of the secondary in massive binaries. 
For an independent mass distribution of the secondary, the 
power law index must be fiatter than 0.3. In particular, an 
independent mass distribution of a Kroupa or Salpeter type 
is strongly excluded. Assuming that the masses of compo- 
nents in a binary are not independent, our results are con- 
sistent with the flat mass ratio distribution. 

The fraction of compact objects. X-ray active in 



LMXBs, on the contrary, is small, fx ~ 10~ , demonstrat- 
ing that LMXBs are extremely rare objects. This result is 
in line with the conclusions of the binary evolution theory 
(e.g. Verbunt & van den Heuvel 1995). 
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APPENDIX A: CATALOGUE OF HMXB IN 
NEARBY STAR-FORMING GALAXIES 

We present a catalogue of compact X-ray sources detected 
within the D25 ellipse for galaxies from our primary sam- 
ple. The full catalogue data includes 1057 sources and is 
currently available at the link shown in the note to the 
Table lAll The latter includes only a small fraction of the 
catalogue. Table \KT\ lists the following information for each 
compact sources: sequence number (col. 2), Chandra X-ray 
Observatory (CXO) source name (JHHMMSS.s-DDMMSS) 
(col. 3), right ascension and declination (J2000.0) in degrees 
(cols. 4 and 5), distance from the galaxy center in arcsec 
(col. 6), number of source counts after background subtrac- 
tion and its statistical error computed as described in Sec. 
3.2 (cols. 7 and 8), decimal logarithm of X-ray flux and lu- 
minosity, measured in the 0.5 — 8keV band as described in 
Sec. 3.2 (cols. 9 and 10), location of the source within the 
-D25 (col. 11). The latter indicates whether the point source 
has been detected in the inner galactic region dominated by 
HMXB population and selected as described in Sect. U (flag 
1), in the outer region dominated by CXB sources (flag 2) 
or in the bulge region dominated by LMXB population (flag 
3). 

This paper has been typeset from a l^jiX/ IM^TjX file prepared 
by the author. 
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Table Al. List of compact sources detected within the D2b eUipse in all galaxies of the resolved sample. 



Galaxy 


ID 


CXO name 


"J2000 


'5.12000 


d 


S 


<7S 


log Fx 


log Lx 


location 








(deg) 


(deg) 


(arcsec) 


(cts) 


(cts) 


(erg cm~^ ^~^) 


(ergs^-^) 


flag 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


NGC 0278 


1 


CXOU J005204. 3+473304 


13.01831 


47.55115 


2.5 


41.7 


9.6 


-14.42 


37.80 


3 




2 


CXOU J005204.3+473258 


13.01823 


47.54965 


3.1 


28.7 


8.4 


-14.58 


37.64 


3 




3 


CXOU J005204.6+473305 


13.01957 


47.55146 


5.2 


50.1 


10.0 


-14.34 


37.88 


3 




4 


CXOU J005203.8+473309 


13.01596 


47.55265 


9.1 


62.5 


10.3 


-14.20 


38.02 


1 




5 


CXOU J005203.4+473254 


13.01429 


47.54855 


11.4 


74.9 


11.0 


-14.15 


38.07 


1 




6 


CXOU J005204.4+473314 


13.01835 


47.55404 


12.8 


30.5 


8.0 


-14.54 


37.68 


1 




7 


CXOU J005203.3+473312 


13.01410 


47.55333 


13.9 


25.1 


7.4 


-14.62 


37.60 


1 




8 


CXOU J005205.8+473302 


13.02457 


47.55056 


16.1 


31.5 


7.8 


-14.54 


37.68 


1 




9 


CXOU J005202.7+473307 


13.01137 


47.55217 


17.1 


1035.4 


39.1 


-13.00 


39.22 


1 




10 


CXOU J005205.5+473249 


13.02299 


47.54716 


17.2 


11.5 


6.7 


-14.98 


37.25 


1 




11 


CXOU J005202.1+473316 


13.00896 


47.55446 


26.1 


16.3 


6.3 


-14.79 


37.43 


1 




12 


CXOU J005205.9+473321 


13.02485 


47.55610 


26.2 


17.9 


6.5 


-14.78 


37.44 


1 




13 


CXOU J005206.0+473347 


13.02540 


47.56320 


49.2 


114.7 


13.5 


-13.98 


38.24 


1 


NGC 0520 


1 


CXO J012435. 1+034731 


21.14660 


3.79199 


2.3 


1100.1 


37.2 


-12.79 


40.17 


1 




2 


CXO J012434.8+034729 


21.14552 


3.79162 


3.6 


50.6 


11.6 


-14.13 


38.84 


1 




3 


CXO J012435. 6+034729 


21.14865 


3.79150 


9.7 


58.1 


10.0 


-14.07 


38.90 


1 




4 


CXO J012434.1+034740 


21.14237 


3.79447 


15.4 


28.0 


7.3 


-14.34 


38.62 


1 




5 


CXO J012433.5+034733 


21.13973 


3.79259 


23.0 


15.0 


5.7 


-14.63 


38.34 


1 




6 


CXO J012433.5+034748 


21.13958 


3.79681 


28.3 


465.8 


24.5 


-13.14 


39.82 


1 




7 


CXOU J012436.5+034659 


21.15243 


3.78320 


40.2 


9.1 


4.7 


-14.88 


38.09 


1 




8 


CXO J012437.6+034715 


21.15702 


3.78772 


42.6 


58.8 


9.8 


-14.07 


38.90 


1 




9 


CXO J012432.7+034806 


21.13660 


3.80198 


48.5 


24.1 


6.8 


-14.45 


38.52 


1 




10 


CXO J012431. 9+034703 


21.13308 


3.78437 


55.1 


13.5 


5.4 


-14.70 


38.26 


1 




11 


CXO J012438.7+034742 


21.16147 


3.79512 


56.0 


16.4 


5.8 


-14.62 


38.35 


1 




12 


CXO J012438.4+034812 


21.16015 


3.80340 


64.0 


13.0 


5.5 


-14.70 


38.27 


2 




13 


CXO J012437.8+034625 


21.15756 


3.77382 


78.6 


45.9 


9.1 


-14.18 


38.79 


1 




14 


CXO J012435.5+034604 


21.14791 


3.76794 


88.3 


20.4 


6.5 


-14.53 


38.44 


2 




15 


CXOU J012441. 5+034607 


21.17299 


3.76872 


128.8 


12.9 


5.7 


-14.70 


38.27 


2 


NGC 1313 


1 


CXO J031818. 2-663004 


49.57604 


-66.50116 


16.8 


385.1 


22.5 


-12.86 


38.45 


1 




2 


CXO J031818.8-663001 


49.57878 


-66.50042 


18.8 


304.0 


20.8 


-12.96 


38.35 


1 




3 


CXO J031820.0-662910 


49.58346 


-66.48642 


48.8 


5137.4 


79.6 


-11.68 


39.62 


1 




4 


CXOU J031821. 2-662858 


49.58867 


-66.48286 


63.6 


15.0 


5.7 


-14.26 


37.05 


1 




5 


CXO J031805.4-663014 


49.52287 


-66.50414 


66.6 


161.2 


14.9 


-13.22 


38.08 


1 




6 


CXO J031806.3-663038 


49.52662 


-66.51073 


73.2 


64.2 


10.4 


-13.63 


37.68 


1 




7 


CXO J031823.7-662834 


49.59902 


-66.47625 


91.7 


19.6 


6.4 


-14.14 


37.17 


1 




8 


CXOU J031829.5-662841 


49.62308 


-66.47809 


108.6 


13.0 


5.2 


-14.32 


36.99 


1 




9 


CXO J031830.5-663135 


49.62741 


-66.52659 


134.0 


25.6 


6.6 


-14.04 


37.27 


2 




10 


CXO J031839.5-662920 


49.66462 


-66.48898 


144.3 


14.2 


5.6 


-14.29 


37.02 


2 




11 


CXO J031818.8-663230 


49.57863 


-66.54172 


157.4 


85.7 


11.1 


-13.51 


37.80 


2 




12 


CXO J031748.7-663043 


49.45354 


-66.51204 


170.1 


32.1 


7.6 


-13.91 


37.39 


2 




13 


CXO J031757.7-663224 


49.49069 


-66.54032 


186.8 


33.1 


7.6 


-13.91 


37.40 


2 




14 


CXO J031747.3-663121 


49.44733 


-66.52265 


192.7 


320.1 


21.3 


-12.92 


38.39 


2 




15 


CXO J031742.5-663123 


49.42745 


-66.52330 


219.4 


132.2 


14.0 


-13.28 


38.03 


2 




16 


CXOU J031735.6-663059 


49.39851 


-66.51654 


250.4 


13.1 


5.4 


-14.28 


37.03 


2 



Note. This table is presented in its entirety in the electronic version (http://www.mpa-garching.mpg.de/~mineo). An abbreviated version 
of the table is shown here for guidance as to its form and content. Full catalogue data includes 11 columns of information for 1057 sources. 
Meanings and units for all columns have been summarized in detail in Appendix A. 



